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In  this  dissertation  we  have  developed  a "transfer-function"  type 
theory  for  linear  time-varying  discrete-time  systems.  Using  this 
framework  in  the  first  part  of  the  dissertation,  we  have  been  able  to 
generalize  much  of  the  existing  polynomial  model  theory.  Specifically 
we  have  treated  polynomial  realization  theory  (FUHRMANN),  polynomial 
factorization  theory,  and  applications  to  feedback  control.  In  the 
second  half  of  the  dissertation,  we  have  treated  the  problems  of 
stabilization,  of  existence  of  stable-proper  factorizations,  and  have 
taken  a cursory  look  at  the  tracking  problem  for  time-varying  systems. 
One  of  our  most  significant  results  is  the  equivalence  of  dynamic  and 
memoryless  state  feedback  as  far  as  the  problem  of  stabilization  is 


concerned . 


CHAPTER  ONE 
INTRODUCTION 


This  dissertation  is  concerned  with  the  study  of  linear 
discrete-time  time-varying  systems.  Time-varying  systems  arise 
frequently  in  practical  applications.  For  example,  time  variation 
could  result  from  change  of  the  mass  and  center  of  gravity  of  an 
aircraft  due  to  fuel-burn,  aging  or  slag  buildup  in  chemical 
reactors,  linearization  of  non-linear  systems  about  a time-varying 
nominal  trajectory,  etc.  The  broad  objective  of  this  dissertation  is 
to  develop  a systematic  theory  for  the  analysis  of  linear  time- 
varying  systems  based  on  matrix-f raction  representations  and  also  to 
apply  this  theory  to  feedback  control  system  design  problems. 

The  earliest  approaches  to  studying  feedback  control  problems 
for  linear  time-varying  systems  were  based  on  linear  quadratic 
optimal  control  theory  and  pioneered  by  Kalman.  In  particular  KALMAN 
[1960]  showed  that  a uniformly  reachable  continuous-time  time-varying 
system  can  be  stabilized  by  a state  feedback  control  law  of  the  form 
u(t)  = -L(t)  x(t).  Here  u(t)  is  the  input  to  the  system  and  x(t)  is 
the  state  of  the  system.  The  gain  matrix  L(t)  may  be  computed  by 
solving  a time-varying  Riccati  differential  equation.  The  reader  may 
consult  the  books  by  JAZWINSKI  [1970]  and  KWAKERNAAK  and  SIVAN  [1972] 
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for  further  results  on  optimal  control/filtering  of  linear  time- 
varying  systems. 

CHENG  [1979]  has  obtained  different  stabilizing  memoryless 
feedback  laws  (based  on  controllability  grammians)  as  have  KWON  and 
PEARSON  [1977,  1978]  (based  on  receding  horizon  optimal  control). 
Both  Cheng  and  Kwon  and  Pearson  work  with  reachable  time-varying 
systems.  Recently,  ANDERSON  and  MOORE  [1981a]  have  been  able  to 
define  the  weaker  notion  of  stabilizability  and  have  obtained 
memoryless  stabilizing  feedback  laws  for  stabilizable  time-varying 
systems.  Their  work  is  based  on  an  optimal  control /fi Iteri ng 

approach  and  specified  in  terms  of  Riccati  difference  equations. 

In  contrast  to  the  above-mentioned  work  which  focusses  on 
obtaining  stabilizing  feedback  laws  based  on  grammi ans/Ri ccati 
equations,  this  dissertation  is  an  attempt  to  develop  a more 
algebraic  theory  for  the  study  of  linear  time-varying  systems.  An 
algebraic  theory  would  provide  more  insight  in  the  study  of  feedback 
control  problems  and  would,  perhaps,  yield  control  laws  that  are 
easier  to  compute.  In  the  past,  there  has  been  some  effort  directed 
at  obtaining  an  algebraic/geometric  theory  for  time-varying  systems 
(see  for  example  WOLOVICH  [1968],  MORSE  and  SILVERMAN  [1972],  KAMEN 
and  HAFEZ  [1979]).  This  work,  however,  treats  very  restricted 
classes  of  linear  time-varying  systems  such  as  i ndex-i nvariant 
systems  or  cyclizable  systems.  This  dissertation  attempts  the  study 
of  linear  time-varying  systems  in  complete  generality. 
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In  Chapters  Two  and  Three  of  this  dissertation,  we  outline  a 
"transfer-function"  type  theory  for  linear  time-varying  systems. 
Attempts  have  been  made  in  the  past  to  develop  such  a theory  (notably 
the  system  function  of  ZADEH  [1950]),  but  these  have  not  met  with 
much  success.  Our  framework  is  specified  in  terms  of  skew 

(noncommutati ve)  rings  of  polynomials,  formal  power  series,  and 
formal  Laurent  series,  all  with  coefficients  in  the  ring  of  time 
functions.  These  skew  rings  have,  in  previous  work,  found 
application  to  the  study  of  linear  time-varying  networks  and  systems; 
see  for  example  articles  by  NEWCOMB  [1970],  SALOVAARA  and  BLOMBERG 
[1973],  YLINEN  [1975],  KAMEN  and  HAFEZ  [1979].  The  rudiments  of  the 
transfer-function  approach  we  develop  here  may  be  found  in  an 
unpublished  paper  of  KAMEN  [1974];  however,  a complete  development  of 
this  approach  is  not  attempted  in  that  paper.  KAMEN  and  KHARGONEKAR 
[1982]  have  pursued  this  approach  and  have  developed  much  of  the 
framework  for  this  transfer-function  theory  for  linear  time-varying 
systems.  Indeed,  this  dissertation  is  a natural  extension  of  their 
work . 

During  the  past  decade,  significant  progress  has  been  made  in 
the  study  of  both  linear  time-invariant  systems  and  systems  over 
commutative  rings  using  polynomial  matrix-fraction  methods  (see  for 
example  the  books  by  ROSENBROCK  [1970],  WOLOVICH  [1974],  KAMEN  and 
ROUCHALEAU  [1984]).  This  approach  has  proven  to  be  useful  in 

tackling  many  system  and  control  theoretic  problems  such  as 

realization,  dynamic  compensation,  regulation  in  the  presence  of 
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disturbances,  etc.  For  details  on  this  work,  the  reader  may  consult 
the  work  of  FUHRMANN  [1976],  ROSENBROCK  and  HAYTON  [1978],  CHENG  and 
PEARSON  [1978],  ANTSAKLIS  [1979],  KHARGONEKAR  [1982],  to  mention  a 
few. 

Given  the  power  and  success  of  polynomial  matrix-fraction 
methods  in  the  linear  time-invariant  system  theory,  it  seems  natural 
to  attempt  to  generalize  this  approach  to  encompass  linear  time- 
varying  systems.  We  do  just  this  in  the  first  half  of  the 
dissertation.  This  generalization  to  time-varying  systems  is 
entirely  nontrivial  because  our  framwork  is  specified  in  terms  of  a 
noncommutative  (skew)  ring  structure  to  incorporate  the  time-variance 
of  our  systems.  Consequently,  in  many  instances  we  are  compelled  to 
use  proof  techniques  that  are  novel  and  entirely  different  from  those 
employed  in  the  study  of  linear  time  invariant  systems  over  both 
fields  and  commutative  rings. 

In  particular,  we  obtain  in  Chapter  Four  a "natural"  state-space 
representation  derived  from  a polynomial  matrix-fraction 
representation  of  the  transfer-function  matrix.  This  realization  is 
a time-varying  analog  of  the  FUHRMANN  [1976]  realization  in  the  time- 
invariant  case.  We  then  investigate  the  relationship  between  system- 
theoretic  properties  of  this  realization  and  algebraic  properties  of 
the  associated  polynomial  matrix  representation.  We  also  examine  the 
problem  of  strict  system  equivalence  and  derive  results  similar  to 
those  obtained  by  FUHRMANN  [1976,  1977]  for  linear  time-invariant 
systems  and  by  KHARGONEKAR  [1982]  for  systems  over  commutative  rings. 
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Many  frequency -doma i n methods  used  to  design  controllers  for 
linear  time-invariant  systems  begin  with  left-  and/or  right-  Bezout 
(coprime)  factorizations  for  the  transfer-function  matrix  of  the 
plant.  As  is  well  known,  any  linear  time-invariant  system  admits 
such  a factorization.  This,  however,  is  not  the  case  for  linear 

time-varying  systems.  In  Chapter  Five,  we  derive  necessary  and 

sufficient  conditions  for  the  existence  of  Bezout  polynomial 
factorizations  for  a linear  time-varying  system.  Moreover,  our 


results  are  construct! ve , and  we  present  a systematic  procedure  for 
obtaining  these  factorizations. 

Following  this,  in  Chapter  Six,  we  use  these  Bezout  polynomial 
factorizations  together  with  the  polynomial  realization  theory  of 
Chapter  Four,  to  study  feedback  control  problems.  In  particular,  we 
derive  an  assignability  result  (which  corresponds  to  being  able  to 
assign  the  closed-loop  system  dynamics)  for  canonical  linear  time- 
varying  systems.  We  illustrate  our  constructive  techniques  by 
designing  a "dead-beat"  controller  for  an  armature-control  led  dc 
motor  with  a time-varying  motor  torque  "constant"  (the  time-variation 
being  due  to  loading  and  heating  effects). 

Chapter  Seven  is  concerned  with  reviewing  some  basic  concepts 
dealing  with  the  stability  of  linear  time-varying  systems,  and  with 
translating  these  concepts  into  our  framework.  In  Chapter  Eight  we 
examine  in  detail  the  problem  of  stabilizing  a linear  time  varying 
system  E.  In  particular  we  introduce  the  notion  of 


asycontrol lability  which  is  equivalent  to  being  able  to  stabilize  Z 
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via  dynami c state  feedback.  ANDERSON  and  MOORE  [1981a]  have  defined 
a notion  of  stabilizability  for  linear  time-varyi ng  systems  which  is 
equivalent  to  being  able  to  stabilize  E via  memoryl ess  state 
feedback.  One  of  the  most  significant  results  in  this  dissertation 
is  the  equivalence  of  stabilizability  and  asycontrollablity . This 
result,  in  particular,  implies  that  dynamics  in  state  feedback  buy 
nothing  extra  as  far  as  the  problem  of  stabilization  is  concerned. 

Of  late  the  use  of  stable-proper  factorizations  in  the  analysis 
and  design  of  linear  time-invariant  control  systems  has  become 
increasingly  popular,  one  advantage  being  that  properness  of 
controllers  is  automatic.  See  for  example  the  work  of  DESOER  et  al . 
[1980],  VIDYASAGAR  [1978],  and  SAEKS  and  MURRAY  [1981].  In 

Chapter  Nine  we  examine  in  detail  stable-proper  factorizations  for 
time-varying  systems.  Following  this,  in  Chapter  Ten  we  investigate 
the  role  of  these  factorizations  in  feedback  control  problems.  In 
particular  we  formulate  the  problem  of  Tracking  with  Internal 
Stability  (TIS)  for  time-varying  systems,  and  show  using  stable- 
proper  factorizations  that  the  TIS  problem  can  be  solved  if  and  only 
if  a particular  linear  matrix  equation  over  a skew  ring  admits  a 
solution. 

Finally,  in  Chapter  Eleven,  we  make  some  concluding  remarks  and 
discuss  some  open  problems  in  the  area  of  linear  time-varying 


systems . 


CHAPTER  TWO 

PRELIMINARY  DEFINITIONS  AND  RESULTS 


In  this  chapter,  we  first  establish  some  notation  and  state  some 
preliminary  definitions  and  results  on  linear  time-varying  systems. 
We  also  prove  a proposition  on  the  existence  of  a "dead-beat"  control 
law  for  reachable  time-varying  systems. 

As  Z = set  of  integers  and  R = field  of  real  numbers,  let  A 
denote  the  R-1 inear  space  of  all  functions  from  Z into  R.  With  the 
operations  of  pointwise  addition  and  multiplication,  it  is  easy  to 
verify  that  A forms  a commutative  ring  with  identity  1,  where 
l(k)  = 1 for  all  k in  Z.  Of  central  importance  in  this  entire  theory 
is  the  right -shift  operator  a defined  by 

(oa)(k)  = a(k-l),  for  all  k in  Z. 

With  the  shift  operator  a (which  is  a ring  automorphism  on  A),  the 
ring  A is  called  a difference  ring.  A subring  BCA  is  called  a 
difference  subring  of  A if  o(B)  = B. 

(2.1)  EXAMPLE.  Some  examples  of  difference  subrings  are  given 
bel ow. 
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PER(N)  :=  (a  in  A : a is  periodic  with  period  N}. 

oo 

PER  :=  PER(i)  = {set  of  all  periodic  time  functions}. 
i=l 

00 

A (Z)  :=  {set  of  all  bounded  time  functions}. 

R[k]  :=  {set  of  all  time  functions  that  are  evaluations  of 
polynomials  in  time}.  Q 

Let  denote  the  difference  subring  of  A consisting  of  all 
functions  with  support  bounded  on  the  left,  i.e.,  for  any  a in  A+, 
there  exists  an  integer  k such  that  a(k)  = 0 for  k < k . 

(2.2)  DEFINITION.  Let  m and  p be  positive  integers.  An  m-input 
p-output  linear  causal  time-varying  input/output  map  f is  an  R-1 inear 
map 


f 


such  that  if  u(k)  = 0 for  k < k,,,  for  some  u in  then  f(u)(k)  = 0 

foM<_^.  □ 


It  is  well  known  that  for  any  input/output  map  f,  there  exists 
p X m matrix  function  W^(i,j)  such  that  for  any  u in  a’J, 

f(u)(i)  = W^(i,j)u(j). 


a 
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The  matrix  function  is  called  the  unit-pulse  response  function 
associated  with  the  input/output  map  f.  Note  that  by  causality, 

W^(i,j)  is  not  defined  for  i < j. 

Our  next  concept  is  the  notion  of  a system. 

(2.3)  DEFINITION.  Let  B be  a fixed  difference  subring  of  A 
containing  1.  An  m-input  p-output  n -dimensional  linear  time-varying 
system  over  B is  a quadruple  E = (F,G,H,J)  of  matrices  over  B where  F 
isnxn,  Gisnxm,  Hispxn,  and  J is  p x m.  | [ 

For  any  matrix  M over  B,  define  its  conjugate  M’*’  by  M'''(k)  = 
M'(-k)  where  ' denotes  the  transpose.  The  dual  system  I'''  of 

^ = (F,G,H,J)  is  defined  by  E'*’  = (F''',H^,G''’,J''’) . This  time-reversal 

is  an  essential  part  of  the  natural  notion  of  duality  for  time- 

varying  systems. 

With  a system  E = (F,G,H,J),  we  shall  associate  the  dynamical 
equations 


x(j  + 1)  = F(j)x(j)  + G(j)u(j), 
y(j  + 1)  = H(j)x(j)  + J(j)u(j), 

where  x(j),  y(j),  and  u(j)  have  the  usual  interpretation.  In  the 
above  definition  of  a system,  it  is  important  to  observe  that  by 
selecting  the  difference  subring  B,  we  can  restrict  our  attention  to 
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a particular  class  of  systems.  For  example,  we  can  study  the  class 
of  linear  time-varying  systems  with  bounded  coefficients  by  choosing 
B = a”(Z)  = the  set  of  all  bounded  functions  from  Z into  R. 

Let  E = (F,G,H,J)  and  E = (F,G,R,i3)  be  two  m-input,  p-output, 
n-dimensional  systems  over  BCA.  Then,  E and  E are  said  to  be 
(algebraically)  B-i somorphi c if  and  only  if  there  exists  an  n x n 
matrix  T over  B,  with  inverse  T"^  over  B,  such  that 


As  is  well  known,  two  isomorphic  systems  E and  E are  related  via  a 
coordinate  transformation  x(k)  = T(k)x(k)  of  their  states. 

Also,  the  unit  pulse  response  function  Wj.  associated  with  the 
system  E = (F,G,H,J)  is  given  by 


The  input/output  behaviour  of  E is  described  by  its  input/output  map 
fj;  where 


F = FT,  R = HT 


g = (a"^T‘^)G,  J = J 


Wj.(i,j)  = J(j) 


H(i)F(i-l)F(i-2)  ...  F(j+l)G(j),i>j 


not  defined 


.i<j 


1^2;(u)(i)  = Wj.(i,j)u(j) 
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Given  an  input  map  f : a'J  ^ aJ,  a system  E = (F,G,H,J)  is  said 
to  be  a real i zati on  of  f if  and  only  if  f = fj..  For  further  results 
on  realizability,  we  refer  the  reader  to  WEISS  [1972],  EVANS  [1972], 
and  FERRER  and  KAMEN  [1984]. 

A system  E = (F,G,H,J)  or  the  pair  (F,G)  over  A is  said  to  be 
reachable  in  N steps  if  and  only  if  there  exists  a positive  integer  N 
such  that  for  any  j in  Z and  any  x in  R*^,  there  exists  an  input 
sequence  u(j-N),  u(j-N+l),  ...,  u(j-l)  which  drives  E from  x(j-N)  = 0 
to  x(j)  = X.  The  dual  notion  of  observability  in  N steps  has  the 
obvious  system-theoretic  interpretation. 

Also,  a system  E is  observable  in  N steps  if  and  only  if  its 
dual  E+  is  reachable  in  N steps.  A system  E is  said  to  be  canonical 
if  and  only  if  it  is  both  reachable  and  observable  in  N steps. 

Let  R^  denote  the  N-step  reachability  matrix 

Rn  :=  [G  F (aG)  ...  F(oF)  ...  ( a^'^F) ( a^'^G) ] . 

WEISS  [1972]  has  obtained  the  following  characterization  of 

reachability. 

(2.4)  LEMMA.  The  pair  (F,G)  over  A is  reachable  in  N steps  at  all 

times  if  and  only  if 


rank  Rm(j)  = n 


for  all  j in  Z, 
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i .e. , 1f  and  only  if  R|^  is  right-invertible  over  A.  [ [ 

(2.5)  REMARK.  For  time-varying  discrete-time  systems,  it  can  happen 
that  a pair  (F,G)  is  reachable  in  N > n steps  (see  example  (2.6)). 
This  is  due  to  the  lack  of  a "Cayley-Hamilton"  type  theorem  in  this 
setting.  Q 

(2.6)  EXAMPLE.  Consider  the  pair  (F,G)  over  A where  F(k)  = 1 for 
all  k in  Z,  and  G(k)  = 1 if  k is  even  and  G(k)  = 0 otherwise.  The 
pair  (F,G)  is  easily  seen  to  be  reachable  in  two  steps  but  not  in  one 
step.  Q 


In  some  instances,  one  may  be  interested  in  a slightly  different 
notion  of  reachability.  Let  T.  = (F,G,H,J)  be  a system  over  a 

difference  subring  B of  A.  Then,  E is  said  to  lDe  B-reachable  in  N 
steps  if  and  only  if  Rj^j  is  right-invertible  over  B.  For  example,  if 
B = Ji“(Z),  this  notion  of  reachability  is  equivalent  to  requiring 
uniform  boundedness  (with  respect  to  j)  of  the  inputs  u(j-N), 
U(j-N+1),  ...,  in  the  definition  of  reachability  given  earlier.  In 
fact  we  have  the  following  result. 

(2.7)  PROPOSITION.  Let  I = (F,G,H)  be  a time-varying  system  over 
£~(Z).  Then,  E is  £“(Z)-reachable  if  and  only  if  there  exists  a real 


number  e such  that 
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> e > 0 


PROOF.  Suppose  D is  £ (Z)  reachable.  Then,  R|^  is  right-invertible 
with  inverse  U over  £"(Z).  Since  U is  over  £"{Z),  this  implies  that 
there  exists  a real  number  M such  that  all  the  nxn  minors  u^-  of  U are 
bounded  by  M.  Let  r^  denote  the  nxn  minors  of  R[^.  Then,  by  the 
Cauchy-Binet  formula  (see  GANTMACHER  [1959])  it  follows  that  for  any 
time  k 


The  last  inequality  in  particular  implies  that  there  exists  an  e > 0 


Conversely,  suppose  det(R^RJ^) (k ) > e > 0 for  all  time  k. 

This  implies  that  is  right-invertible  with  right  inverse 


det(R^R')(k)  = Z r^(k) 


det(R^U)(k)  = 1 = S r.(k)u.(k)  < M E |r.(k)| 


such  that  Z r:(k)  > e which  in  turn  implies  that  det(R|^Rj^) (k ) > e. 


It  is  now  clear  that  V must  be  bounded 


since  det(Rj^Rj^)  > e. 

This  completes  the  proof. 


□ 


We  shall  also  need  the  following  result  on  reachability. 
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(2.8)  PROPOSITION.  Let  I = (F,G,H,J)  over  A be  reachabl  e in  N steps 
at  all  times.  Then,  there  exists  an  n x m matrix  L over  A such  that 
(F-GL)  1s  g-ni1 potent,  1.e.,  for  some  integer  q, 

(F  - GL)a(F  - GL)  ...  a^(F  - GL)  = 0. 

PROOF.  Consider  the  linear  time-varying  system  described  by 

x(k  + 1)  = F(k)x(k)  + G(k)u(k)  . 

Since  the  pair  (F,G)  is  reachable  in  N steps,  it  is  also  controllable 
to  the  origin  in  N steps,  i.e.,  any  initial  state  x(0)  in  can  be 
driven  to  zero  final  state  in  N steps.  Therefore,  there  exist 
controls  u(0),  u(l),  ...,  u(N  - 1)  in  such  that 

/ 

x(l)  = F(0)x(0)  + G(0)u(0) 
x(2)  = F(l)x(l)  + G(l)u(l) 

x(N)  = F(N  - l)x(N  - 1)  + G(N  - l)u(N  - 1)  = 0 . 

Let  C be  the  subspace  consisting  of  all  initial  states  x(0)  in 
R^  that  can  be  driven  to  zero  final  state  in  i -steps  or  less. 
Cl  early 


V^C  V^C  ...iZ  = r""  . 
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We  now  select  a linearly  independent  set  of  vectors  bp  b2,  ....  b^ 
that  span  Vp  We  then  extend  this  set  to  form  a basis  for  V^.  In 
this  manner,  we  find  a basis  B = {bp  b^}  for  = R^.  Suppose 

some  b^  in  B is  in  but  not  in  V^_p  For  the  initial  state 
x(0)  = bp  we  can  therefore  find  a control  sequence  u^O),  u''(l), 
u^(N  - 1)  where  u’(k)  = 0 for  k > t which  drives  x(0)  to  the 
final  state  x(N)  = 0.  In  this  manner,  we  determine 


{u’(k)  ; k = 0,  1,  ...,  N - 1 ; i = 1,  2 n}  . 


Let  B be  the  n X n invertible  matrix  over  R whose  i-th  column  is 
bi . Let  Uq,  Up  ...,  U|y|_2  be  m X n matrices  over  R where  the  i-th 
column  of  U|^  is  u^(k).  Also,  recursively  define  the  matrices  X^,  X^, 
. . . , X|^  by 


^k+1  “ \ 

It  then  follows  that  Xf^  = 0. 

Suppose  5 is  in  Ker  (X^-).  This  means  that  the  initial  state 
x(0)  = B C can  be  driven  to  zero  final  state  in  i steps.  Therefore, 
B 5 is  in  Vp  and  then  by  our  judicious  choice  of  controls,  it 


follows  that 
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U|^  5 = 0 , for  k > i 

Specifically,  5 = 0 or  Ker(X^)  _ Ker(U^.  )•  Hence  we  can  solve  the 
N linear  equations 

U.  = -L.X.  , i = 0,  1,  ....  N - 1 

for  L^- . Observe  now  that 

\+l  = \ 

= (F(k)  - G(k)  L^)  X^  . 


It  then  follows  that 

0 = X,^  = (F(N-l)  - G(N-l)  ...  (F(0)  - G(O)Lo)  X^  . 

However,  Xq  = B is  invertible.  Consequently, 

(F(N-l)  - G(N-l)  L,^_i)(F(N-2)  - G(N-2)  Lf^_2) . . . (F(0)  - G(0)  L^)  = 0 . 

Since  (F,  G)  is  reachable  in  N steps  at  all  times,  we  can  repeat 
the  above  argument  and  find  matrices  ...  , L2i^_i  such  that 


(F(2N-1)  - G(2N-1)  L2n_i)  ...  (F(N)  - G(N)  L^)  = 0 . 
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In  this  fashion,  we  find  matrices  L|^  for  each  integer  k.  Define  the 
m X n matrix  L over  A by  L(k)  :=  L|^.  It  is  then  apparent  that 

(F  - GL)a(F-GL)  ...  (F-GL)  = 0 , 

i.e.,  (F  - GL)  is  a-nilpotent.  EH 

(2.9)  REMARK.  It  is  interesting  to  note  that  the  above  proposition 
implies  the  existence  of  a "dead-beat"  control  law  for  linear  time- 
varying  systems  that  are  reachable  in  N steps.  More  precisely,  with 
the  matrix  L as  described  in  Proposition  (2.8),  it  follows  that 
starting  from  any  initial  state  ? in  at  any  initial  time  k^,  the 
state  trajectory  x(k)  of 

x(k+l)  = (F  - GL)(k)  x(k) 


is  zero  for  k > k^  + q. 


□ 


CHAPTER  THREE 

THE  TRANSFER-FUNCTION  FRAMEWORK 


In  this  chapter,  we  describe  in  detail  the  elements  of  a 
"transfer-functi on"  type  approach  to  linear  time-varying  systems 
introduced  by  KAMEN  and  KHARGONEKAR  [1982].  Much  of  the  remainder  of 
this  dissertation  is  based  on  this  theory. 

The  commutative  rings  of  polynomials,  power  series  and  formal 
Laurent  series  all  with  coefficients  in  the  reals  R,  play  a central 
role  in  the  transfer-function  theory  of  linear  time-invariant 
systems.  For  time-varyi ng  systems  the  analogous  objects  are  skew 
(noncommutati ve)  rings  with  coefficients  in  the  ring  of  time 
functions. 

More  precisely,  with  z equal  to  an  indeterminate,  let  A((z"^)) 
denote  the  set  of  all  formal  Laurent  series  of  the  form 

rl-N  ^"''“r  ’ “r  in  A 

Note  that  the  coefficients  above  are  written  on  the  right.  This  is 
because  we  will  now  impose  a non-commutati ve  ring  structure  on 
A((z"^)).  With  the  usual  addition,  and  with  multiplication  defined 
by 
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(3.1) 

az  = z(aa)  , a in  A , 

where  (aa)(k)  = a(k-l),  A((z“^))  is  a noncommutati ve  ring  with 

identity,  called  the  skew  ring  of  formal  Laurent  series  over  A. 
There  are  two  important  subrings  of  A((z"^)):  The  skew  ring  of 

polynomials  A[z]  and  the  skew  ring  of  formal  power  series  A[[z”^]]. 
These  have  the  obvious  definitions. 

The  noncommutati ve  multiplication  in  A((z'^))  defined  above 
captures  in  a very  natural  way  the  time-variance  of  our  systems,  and 
thus  plays  a central  role  in  this  entire  theory.  This  is  illustrated 
by  the  following  example  and  will  become  more  evident  in  the 
remainder  of  this  chapter. 

(3.1)  EXAMPLE.  Consider  the  following  single-input  single-output 
time-varying  difference  equation: 

a(k+l)y(k+l)  = u(k) 

The  indeterminate  z~^  will,  as  in  the  time-invariant  theory, 
represent  a delay  operator.  The  above  equation  can  be  written  (we 
show  this  more  formally  later)  as 
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z(ay)  = u 

We  can  also  write  the  above  difference  equation  as 

where  3 = (o  ^a) , and  this  can  be  represented  in  the  "frequency- 
domain"  by  zy  = ju,  or,  multiplying  on  the  left  by  3,  3zy  = u. 
Comparing  this  with  the  previous  representation,  we  obtain  z^^  = 3z  = 
(a"^a)z  which  is  precisely  the  manner  in  which  we  have  defined  our 
non-commutati ve  multiplication  in  A((z“^)).  ED 

Define  a projection  map 

TT  : A((z'^))  > A((z‘^))  ^ z'^^  . 

For  any  a in  A((z~^)),  let  (a)^  :=  a-ir(a)  = the  polynomial  part 
of  a.  By  (a)^,  we  shall  mean  the  constant  coefficient  of  a,  and,  a 
is  said  to  be  strictly  proper  if  and  only  if  Tr(a)  = a. 

Given  an  r x r skew  polynomial  matrix  (i.e.,  a matri'x  with 
entries  in  A[z]), 

i 

0 = 2 z'o..  , 

i=0  ^ 
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the  degree  of  Q,  written  deg(Q),  is  the  largest  integer  q such  that 
Oq  * 0.  Further,  Q is  said  to  be  monic  if  Qq  = I (the  r x r identity 
matrix),  and,  Q is  said  to  be  ri ght-i nverti b1  e if  and  only  if  there 
exists  an  r X r matrix  over  A((z"^))  such  that  = I.  We  shall 
need  the  following  result  on  invertibil ity. 

(3.2)  PROPOSITION.  Let  Q be  an  r x r polynomial  matrix.  Then,  Q is 
right-i nverti bl e if  and  only  if  there  exists  an  r x r polynomial 
matrix  T such  that  QT  is  monic.  Further,  T can  be  chosen  such  that 
deg(Q)  = deg(QT). 

PROOF.  Suppose  Q is  right-invertible,  i.e.,  there  exists  a in 
ATxr  ^^2  1))  such  that  Qi|;  = I.  Let  deg  Q = d.  We  can  now  write 

z^I  = Q il' 

= Q(^z‘^)^  + Qtt(4'Z^) 

Notice  that  deg  (Qir(i|)z^))  < deg  Q = d.  Therefore  the  highest  degree 
term  of  Q(i|^z^)^  is  z^I.  Choosing  T = (ij^z^)^  which  is  polynomial , 
proves  the  necessity.  Notice  that  deg  QT  = deg  Q = d. 

Now  assume  that  there  exists  a T in  A''^''[z]  such  that  QT  is  a 
nionic  rxr  polynomial  matrix.  We  can  then  do  right  division  of  I by 
QT  and  find  a <p  1n  such  that  Q T <|)  = I,  which  implies 
that  Q is  right-invertible.  Q 
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The  degree  condition  deg(Q)  = deg  (QT)  is  merely  a technical 
fact  which  will  be  useful  in  proving  several  later  results.  An 
analogous  result  holds  for  left-invertibility  of  Q.  In  general, 
1 eft-i nverti bi 1 i ty  is  not  equivalent  to  right-invertibility,  the 
pathology  being  due  to  the  skew  nature  of  our  rings  (see 

Example  (3.3)).  We  shall  almost  always  deal  with  polynomial  matrices 
that  are  both  left-  and  right-  invertible,  in  which  case  we  shall 
call  them  invertible  to  avoid  use  of  cumbersome  prefixes. 

The  following  examples  illustrate  the  skew  multiplication  (3.1) 
in  our  rings  and  contrast  this  with  multiplication  in  R[z]. 

(3.3)  EXAMPLE.  Define  a^,  og,  u in  A by 


o(l-a^),  p :=  OgOj^ 


Also  define  q(z),  hj(z),  h2(z)  by 

q(z)  = a^z  + a^,  h^(z)  = a^z  + a^(a^)  , h^{z)  = pz  - a(y)  . 


Using  multiplication  in  the  skew-ring  A[z]  defined  by  (3.1)  it  is  an 
easy  computation  to  show  that 


q(z)h^(z)  = z , h2(z)q(z)  = 0 
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Thus,  q(z)  is  right-invertible  but  left-invertible, 
that  does  not  occur  in  R[z], 


a phenomenon 

□ 


(3.4)  EXAMPLE.  Define  3 in  A by 


3(k) 


k is  even 

k i s odd 


Consider  the  scalar  poynomial  q(z)  :=  3z  + 1.  It  is  easy  to  verify 
that 


q(z)h(z)  = h(z)q(z)  = 1 

where  h(z)  = -3z+l.  Thus,  q(z)  is  invertible  and  its  inverse 
q ^(z)  = h(z)  is  polynomi al . This  situation  also  is  peculiar  to  the 
skew-ring  A[z].  []]] 

We  now  describe  a transfer-function  approach  to  time-varying 
systems  based  on  the  skew  rings  defined  earlier.  All  proofs  are 

omitted;  they  can  be  found  in  KAMEN,  KHARGONEKAR  and  POOLLA  [in 
press]. 

Again,  let  A+  denote  the  subring  of  A consisting  of  all 
functions  a : Z R with  support  bounded  on  the  left.  Let  A denote 
the  unit  pulse  at  the  origin,  i.e.. 
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if  k=0 
otherwi se 


Given  any  u in  A^,  the  (generalized)  z-transform  of  u written  U(z)  is 
defined  to  be  the  skew  Laurent  series 

(3.5)  U(z)=  Sz“^u(r)A 

This  generalized  z-transform  is  simply  the  usual  z-transform 
(imbedding  R((z"^))  in  A((z'^)))  multiplied  by  A. 

Lef  f be  an  input/output  map,  and  let  denote  the  unit-pulse 
response  function  associated  with  f.  For  each  integer  r > 0,  define 
a p X m matrix  over  A by 

Wp(k)  = W^(r  + k,  k)  , k in  Z 


(3.6)  DEFINITION.  The  (formal)  transfer-function  matrix  W^(z) 
associated  with  the  input/output  map  f is  the  p x m matrix  over 
A[[z~^]]  defined  by 


Wf(z) 


z 


□ 
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This  definition  of  transfer-function  for  time  varying  systems 
will  be  seen  by  the  following  Propositions  to  be  the  natural 
definition  to  capture  the  time-varying  behaviour  of  our  systems. 

(3.7)  PROPOSITION.  Let  f be  an  input/output  map.  Let  y be  the 

output  resulting  from  the  input  u in  a'^.  Let  Y(z)  and  U(z)  denote 

the  (generalized)  z-transforms  of  y and  u respectively.  Then 

(3.8)  Y(z)  = W^(z)U(z)  . □ 

Note  the  close  resemblance  of  (3.8)  to  the  time-i nvari ant 
transfer-function  theory.  Proposition  (3.7)  is  a result  one  would 
desire  of  any  definition  for  transfer-function . This  analogy  to 

time-invariant  systems  is  further  illustrated  by  the 

(3.9)  PROPOSITION.  Let  I = (F,G,H,J)  be  a linear  time -varying 

input/output  map  f^.  Then,  the  transfer-function  matrix 

A 

Wj.  associated  with  f^.  is  given  by 

(3.10)  Wj,(z)  = H(zl  - F)“^  + J . □ 

Despite  the  close  resemblance  these  two  results  bear  to  the  time- 
invariant  theory,  it  must  be  emphasized  that  (3.8)  and  (3.10)  are 
computed  via  the  skew  (noncommutati ve)  multiplication  defined 
earlier.  For  instance,  (zI-F)"^  is  determined  by  the  formula 
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(zI-F)'^  = z'^  + (cF)z"^  + (aF)(a^F)z"^  + ...  . 

We  conclude  this  chapter  by  relating  polynomi al  factorizations 
of  transfer-function  matrices  to  collections  of  input/output 
difference  equations.  Consider  the  collection  of  input/output 
difference  equations  with  time-varying  coefficients  described  by 


q r 

i?0  )x(k+i ) = R.  (k+i  )u(k+i ) 

(3.11) 

t 

y(k)  = p.(k+i)x(k+i)  . 


Here,  u 
matri ces 
matri ces 


e X e y e A^,  and 

of  appropriate  dimension  over  A. 


the  Q^, 
Defi ne 


, and  R^-  are 
the  polynomial 


Q(Z)  = (io  z'Q,  , R(Z)  . 


r 

i?0 


z’r. 

1 


P(z) 


iio 


z""  P. 
1 


and  assume  that  Q(z)  is  invertible.  Then,  it  is  an  easy  computation 
to  verify  that  the  transfer-function  matrix  associated  with  the 
input/output  map  defined  by  (3.11)  is  given  by 


W^(z)  = P(z)0'^(z)R(z) 


(3.12) 


27 


Conversely,  given  an  input/output  map  f whose  associated  transfer- 

A 

function  matrix  W^(z)  admits  a polynomial  factorization  (as  in 
(3.12)),  one  can  readily  derive  a collection  of  input/output 
difference  equations  (as  in  (3.11))  that  correspond  in  a natural  way 
to  the  particular  factorization  (3.12). 

For  linear  time-invariant  systems  this  correspondence  was  first 
observed  by  ROSENBROCK  [1970].  Polynomial  factor! zat  ions  of 

transfer-function  matrices  are  investigated  extensively  in  the  next 
three  chapters. 
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CHAPTER  FOUR 

POLYNOMIAL  REALIZATION  THEORY 

We  will  now  consider  polynomial  matrix  fraction  representations 
of  a given  transfer  function.  Let  P,  Q,  R,  and  S be  polynomial 
matrices  sucn  that  Q is  invertible.  Let  f be  an  input/output  map 
with  the  associated  transfer  matrix 

W^(z)  = PQ'^R  + S 

For  time-invariant  systems  over  fields,  FUHRMANN  [1976]  gave  a 
realization  for  f in  terms  of  the  polynomial  matrices  P,  0, 
R,  and  S.  For  time-invariant  systems  over  arbitrary  commutative 
rings,  KHARGONEKAR  [1982]  has  obtained  corresponding  results.  We  now 
proceed  to  derive  results  for  time -vary i nq  systems  analogous  to  those 
of  Fuhrmann  and  of  Khargonekar.  We  first  develop  the  machinery  with 
which  we  can  obtain  these  natural  realizations  for  time-varying 
systems . 

Let  Q be  an  rxr  (skew)  polynomial  matrix,  invertible  over 
A((z"^)).  Define  a right  A-module 

Xq  :=  {x  in  A*'[z]  : Q~^x  is  strictly  proper} 
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Define  a right  A-1 inear  projection  map 

: A>'[z]  Xg  : X ->■  Qtt(O'^x) 

where  tt(Q"^x)  is  the  strictly  proper  part  of  Q'^x.  Clearly,  irg  is 
surjective.  The  map  wg  corresponds  to  viewing  polynomials  in  A’'[z] 
modulo  Q. 


(4.1)  PROPOSITION.  Xg  is  a finitely  generated,  free  right  A-module 

(j.e.,  Xg  is  isomorphic  as  a right  A-module)  to  A*^  for  some  i nteger 

n . 


PROOF.  The  proof  of  this  technical 
therefore  put  it  in  Appendix  A. 


fact  is  rather  long,  and  we  have 


□ 


(4.2)  EXAMPLE.  Let  us  consider  the  monic  rxr  polynomial  matri; 


Q(z)  = z"  + + ...  + zQ^  + Og  . 


Let  denote  the  j-th  column  of  the  rxr  identity  matrix.  Notice 
that  for  i = 0,  1,  ...,  n-1;  j = 1,  2,  ...,  r,  z^^j  is  in  Xg  since 
is  strictly  proper.  In  fact  the  set 


^ ^ ,®.j  • ^ “0>  1»  •••  , n-1  ;j  = 1,  2,  ...  , r} 


forms  a basis  for  X^,  and  X = A 


nr 


□ 
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Suppose  Q(z)  is  an  (not  necessarily  monic)  rxr  polynomial  matrix 
of  degree  n.  It  follows  from  the  proof  of  Proposition  (4.1)  that  the 
set 


(4.3) 


{TT^(z^ej.)  : i = 0,  1,  ...,  n-1;  j = 1,  2,  ...,  r} 


where  is  the  j-th  column  of  the  rxr  identity  matrix,  generates 

Xg.  The  difficulty  however,  is  to  extract  a basis  for  Xg  which  we 
shall  require  to  obtain  our  natural  realizations  (see  Remark  (4.5) 
and  Theorem  (4.6)).  The  following  example  illustrates  what  can 
happen  if  0(z)  i s not  moni c . 

(4.4)  EXAMPLE.  Let  q(z)  = 3z^  + z^  where  B(k)  = 1 if  k is  even,  and 
6(k)  = 0 otherwise.  Notice  that  q“l(z)  = z'^  - z-l(a0).  We  compute 
a set  of  generators  for  X^  using  (4.3): 


= q^(q‘^)  = qq’^  = 1 , 


-1 


%(z)  = qi^Cq'^'z)  = qz'^  = 3z^  + z 


^q(z  ) = 0 


Thus,  Xq  is  generated  by  {1,  3z  + z}  which  also  happens  to  be  a 
basis  for  Xq,  and  Xq  = A^.  Q 
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(4.5)  REMARK.  The  right  A-module  Xg  will  serve  as  the  state  space 
for  our  natural  realizations.  Proposition  (4.1)  tells  us  that  for 
any  polynomial  matrix  Q,  our  state  space  Xg  = a"  for  some 
integer  n.  This  means  that  the  dimens i on  of  our  natural  state  space 
(see  Theorem  (4.6))  does  not  vary  with  time:  a comforting  fact.  It 
is  shown  in  the  proof  of  this  proposition  that  the  dimension  of  Xg  is 
given  by  the  rank  of  a certain  matrix  which  can  be  computed  for  any 
given  Q. 

We  can  also  view  Xg  as  a right  A[z]-module  by  defining  right- 
multiplication  by  z in  Xg  as 

X * z = irg(xz)  , for  any  x in  Xg 

This  right  A[z]-module  structure  is  fundamental  for  obtaining  the 
natural  realizations  alluded  to  at  the  beginning  of  this  chapter. 

We  are  now  in  a position  to  state  the  following  result. 

(4.6)  THEOREM.  Let  P,  Q,  R,  and  S be  pxr,  rxr,  rxm,  and  pxm 
polynomial  matrices  such  that  Q is  invertible.  Let  f be  an  input/ 
output  map  with  the  associated  transfer  matrix 

W^(z)  = PQ'^R  + S . 


Define  the  maps  (p,  r,  and  by 
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<}>  : Xp  ; X > TT^(xz) 

r : a"’  -»■  Xq  : u > Hp(Ru) 

4-  : Xq  ^ aP  : X > (PQ"^xz)^  . 

— L»  - . and  H be  the  matrix  representations  of  (|),  r,  and  if) 

respectively  with  reference  to  some  (fixed)  basis  of  Xp.  and  Ipf. 

— •"  i ^ • Then  Z(P,  Q,  _R , $)  :=  (F,  G,  H,  J)  is a 

real ization  of  f. 


PROOF.  Let  {b^,  b2,  b^}  be  a basis  for  Xq,  Any  element  x in  Xq 

has  the  unique  representation  x = E b.  a.,  a.  in  A,  since  Xn 

i=l  1 1 1 0 

is  a free  right  A-moduie  (see  Proposition  4.1).  We  will  write  this 
as 


X = B 


where  B - [b^  b2  ...  b^];  using  this  notation,  the  matrices  F,  G,  and 
H with  respect  to  the  basis  {b^,  b2.  ...,  b^}  are  determined  by 
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BF  = 4,(B)  , 

(4.7)  BG  = r(I  ) , 

m 

H = tj;{B)  , 

where  is  the  mxm  identity  matrix  (over  A). 

Notice  that  the  maps  <t>  and  ij;  and  A-semi  1 i near , i.e., 

<!> 

(xa)  = (j)(x)  a(a) 

'l^(xa)  = i|)(x)a(a) 

for  any  a in  A.  The  map  r is  A-linear.  Let  o denote  composition 
of  maps.  We  now  compute  using  (4.7) 

A^  :=  H(aF)  ...  (a'^'^F)  (0*^3) 

= 'i^(B)(cF)  ...  (a''-^F)(a''G) 

= ip  0 [BF(aF)  ...  (a'^"^F)((j'^"^G)] 

The  last  equality  above  follows  from  the  A-semi  1 i neari ty  of  \|). 
Continuing  our  computation,  we  see  that 
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0 [8F(aF)  ...  (/"^F)(a'^’^G)] 

= ip  0 [<^(B)(aF)  ...  (a^‘^F)(a'^"^G)] 

= (4-  0 <|.)[BF(aF)  ...  (a’^'^F)(/‘^G)] 

where  the  last  equality  was  deduced  from  the  A-semi  1 i nearity  of  <j). 
In  this  manner,  we  see  that 

A,^  = (ip  0 <^'^"^)(BG) 

= 0 0 r)(I^)  . 

Notice  now  that  for  any  a in  A[z],  and  any  x in  A^[z], 

i^Q(^Q(x)a)  = TTQ(xa) 

Therefore, 

A.  = (il'  0 0 r)(i  ) 

^ 171 

= (if-  0 <I>'^"^)(tTq(R))  = (4-  0 <()'^"^)(TrQ(Rz)) 

= (4-)(tTq(Rz'^"^))  = (PQ"^Tr^(Rz^"^)z) 
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= (PQ'^R  . 

From  Proposition  (3.9)  it  follows  that  the  transfer  matrix  Wj,(z)  of 
the  input/output  map  where  Z = (F,  G,  H,  J)  is  given  by 

Wj,  (z)  = H(zl  - F)"^  G + J 
00 

= A,  z"*^  + J 
k=l 

GO 

= z (PO’^R  z*^)^  z'*^  + (PQ-^  R + s) 
k=l  ° 0 

= PQ'^R  + S , 

since  (PQ'^R  + S)  is  proper.  Hence,  Wj,(z)  =W^(z),  i.e. 

Z(P,  0,  R,  S)  = (F,  G,  H,  J)  is  a realization  of  f.  CD 

(4.8)  REMARK.  We  could  have  chosen  a different  basis  {b,,  b^, 

A X C. 

...»  b^}  for  Xg  and  obtained  a different  realization 

A ^ A A 

Z - (F,  G,  H,  J)  for  the  input/output  map  f.  In  this  event,  Z 
and  Z = (F,  G,  H,  J)  are  A-isomorphic.  We  shall  call  any 
realization  of  f obtained  as  in  Theorem  (4.6)  the  Fuhrmann 
real i zati on  associated  with  the  polynomial  matrix  function 
representation  W^(z)  = PQ“^R  + S.  [D 
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We  now  illustrate  Theorem  (4.6)  with  two  examples. 

(4.9)  EXAMPLE.  Let  f be  an  input/output  map  described  by  the 
following  collection  of  input/output  difference  equations 


C(k+2)  + q^C(k+l)  = q^c(k)  = r^u(k+l)  + rgu(k) 
y(k)  = p^5(k+l)  + pQ?(k) 

As  shown  in  Chapter  Three  (see  equations  (3.11)  and  (3.12))  the 

A 

t ransf er-f uncti on  W^(z)  associated  with  f is  given  by  the 

polynomial  representati on 


W^(z)  = PQ-^  = (p^z  + pg)(z^  + zq^  + qQ)'^zr^  + r^^)  . 

where  Pi  = o 

It  follows  from  Example  (4.2)  that  the  set  {1,  z}  forms  a basis 
for  Xq.  Let  us  represent  any  x = oq  + za^  in  Xq  by  the  vector 
[“o  We  first  compute  the  action  of  the  maps  (|),  r,  and  (as 

in  Theorem  (4.6))  on  this  basis  of  Xq: 


<(>(1)  = ■^q(1  • z)  = z 

<}>(z)  = TTj^(z  • z)  = Qtt(Q"^z^)  = z^  - Q(Q"^z^)^ 
= z^  - Q(z)  = -zq^  - qg  , 
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r(l)  = tt^(R  • 1)  = + Tq)  = zr^  + , 

i|»(l)  = (PQ‘h  • z)q  = {(p^z  + Po)(z'^  + ...  )}q  = 

'Kz)  = (PQ'^z-z)q  = {(p^z  + Pq)(1  - z“^a'2(q^)  + ...  )}^ 

_2 

= Pq  - Pi°  (^i) 

Also,  J = (PQ  ^R)q  = p^T]^.  Thus,  the  Fuhrmann  realization  T(P,Q,R)  = 
(F,  G,  H,J)  of  (with  respect  to  the  basis  chosen)  is  given  by 


"0  -qj" 

1 

o 

u 

1 

F = 

» G = 

-’l 

H = [Pq  - Pja”^(q^] 


J = Pir^ 


□ 


(4.10)  EXAMPLE.  Let  f be  an  input/output  map  whose  associated 

A 

transfer  function  W^(z)  admits  the  polynomial  representation  (see 

Example  (4.4)) 


(z) 


Q'^R  = (6z^  + z^) 


+ <'n) 


Observing  that 


<}>(1)  = TT^(1  • z)  = 3 z^  + z 
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4>(3z^  + z)  = ■^g(3z^  + z^)  = 0, 

r(l)  = TT^(zr^  + Pq)  = (ez^  + z)r^  + , 

'i'(l)  = (Q"^*  1 • z)q  = -(cj3)  , 

^(3z^  + z)  = (Q"^(3z^  + z^))q  = 1 , 

we  can  immediately  write  down  the  Fuhrmann  realization  E(I,Q,R)  = 
(F,G,H,J)  of  f as 


1 

o 

1 

-J 

o 

. H = [-{as)  1]  , 

_i  q_ 

-^1- 

□ 

In  the  case  where  0 is  monic,  one  can  readily  obtain  the  general 
form  of  the  Fuhrmann  realization:  Let 


P(z) 

^ i 

= P.z\  Q(z) 

i=0  ^ 

s 

Q.  , R(z)  = E z^R. 

i=0  ^ 

9 

i=0 

and  assume 

that  Oq  = I. 

Let  f 

be  an  input/output 

map  whose 

associ ated 

transfer-functi  on 

ma  t r i X 

A 

W^(z)  admits  the 

polynomi al 

representation 
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W^(z)  = P(z)0“^z)R(z) 


Then,  the  Fuhrman  realization  J:(P,Q,R)  = (F,G,H,J)  of  f is  given  by 


F = 


0 0 


0 

I 


-Qi 


0 ...  I -0^"^ 

q-1 


G = 


O' 

6 


(4.11) 


H = [Pq  P^...P^0...0] 


0 0 . 

0 0 . 


0 I 


a‘^(B2) 

o‘^(B  i) 
-0 


q J 


J = (PQ‘^R)o  . 


where  A“^(z)  = z"^I  + z‘^"^B2  + 83  + ... 


The  Fuhrmann  realization  is  "natural"  in  the  sense  that  control - 
theoretic  properties  of  the  realization  can  be  characterized  in  terms 
of  algebraic  properties  of  t'ne  associated  factorization.  We  have, 
for  instance,  the  following  result  characterizing  reachability  in  N 
steps . 
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(4.12)  THEOREM.  Let  P,  Q,  and  R be  pxr,  rxr,  and  rxm  polynomial 
matrices  over  A such  that  Q 1s  invertible  over  A((z'^)).  Let  f be  an 
input/output  map  with  associated  transfer  matrix 

W^(z)  = PQ'^R  . 

Then,  the  Fuhrmann  realization  £(P,  Q,  R)  = (F,  G,  H)  is  reachable 
in  N steps  for  some  positive  integer  H ir  and  only  if  there  exist 
polynomial  matrices  Y-|  and  such  that 


QYi  + RY2  = I 


PROOF.  Let  i:(P,  Q,  R)  = (F,  G,  H)  be  the  Fuh  rmann  realization  of 

Wf(z)  = PQ  R relative  to  some  (fixed)  basis  {b^,  b2,  b^}  of 

Xq.  We  continue  using  the  notation  used  in  the  proof  of  Theorem 

(4.6):  any  x in  Xq  has  a unique  representation  x = r b.  a. , {a.} 

i=l  1 T T 


in  A,  and  we  shall  write  this  as 


a 


1 


X 


[bj  bj 


Suppose  E - (F,  G,  H)  is  reachable  in  N steps  for  some  integer  N. 

This  means  that  there  exist  matrices  Cq,  C^,  ...,  C(^  ^ 


over  A such 
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that 


rc„  n 


[G|F(0G)|  ...  |F(ctF) 


. (a^-^F)(a^-^G)] 


= I 


Let  ■^q(I)  = BT,  where  T is  an  nxn  matrix  over  A.  It  then  follows 
that 


^q(I)  = B(IT) 


= B [G|F(aG)|  ...  |F(aF)  ...  (a^'“F) (a^'^G)] 


’V 

C^T 


‘V 


= [r(l)|(<t>or)(i)|...|(/-lor)(i)] 


L-W- 


= ^q(RY2) 
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N-1 

where  'i  :=  E z^C.T.  Hence,  it  (I  - RY.)  = 0.  This  implies 
^ i =0  ^ 

that  there  exists  an  rxr  polynomial  matrix  over  A[z]  such  that 

QY^  + RY;,  = 1. 

Now  suppose  that  there  exist  polynomial  matrices  Y^  and  Y2  such 
that  QY^  + RY2  = I.  Then,  for  any  x in  Xq, 


RY2X  = X - QYjx 


Hence,  tTq(RV2x)  = therefore,  write 

B,=  tTq(RY2B)  ' , 

N-1 

where  B is  as  before.  Let  Y2B  = E z^D..  It  then  follows  that 


N-1  . 

B = w (R  E z'd.) 
^ i =1  ^ 


r(i)  + (<(.  0 r)(i)  + ...  + (/-I  0 r)(i) 

B [G|F(aG)|  ...|F...(o^“^F)(a'^"^G)] 


lVu 


Since  {b]^,  b2,  ...,  b^}  is  a basis  for  Xq,  this  implies  that 
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[G|F(aF)|...|F(aF)...(a^-^F)(a^-lG)]  Td^  1 


l'^n-U 


Consequently,  the  system  E = (F,  G,  H)  is  reachable  in  N-steps.  □ 


(4.13)  REMARK.  A dual  version  of  Theorem  (4.12)  characteri zes 

observability  in  N-steps.  We  state  this  without  proof.  The  Fuhrmann 
realization  E(P,  Q,  R)  = (F,  G,  H)  is  observable  in  N steps  if  and 

only  if  there  exist  polynomial  matrices  and  Y/|  such  that 

Y3O  + Y4P  = I . □ 

It  is  a simple  calculation  to  verify  that  the  Fuhrmann 
realization  associated  with  the  polynomial  matrix  representation 
W^(z)  = H(zl  - F)"^G  is  itself  A-isomorphic  to  E = (F,  G,  H). 

Consequently,  from  Theorem  (4.12)  we  immediately  have  the  following 
analog  of  the  familiar  time-invariant  theory  result: 


(4.14)  COROLLARY.  The  pair  (F,  G)  is  reachable  in  N steps  if  and 
only  if  there  exist  polynomial  matrices  Y^  and  Yg  such  that 


(zl  - F)Y^  + GY2  = I 


□ 


44 


Let  Q and  Q be  any  rxr  and  IkZ  invertible  polynoiniai 

Miatricas.  We  ^ill  now  give  necessary  and  sufficient  conditions  for 
Xq  and  Xp  to  be  ACz]-modu1e  isonorphic.  Tor  the  case  A = a fie’  i, 
A(z)  is  the  usual  ring  of  polynomials,  this  issue  was  resolved  by 
FURHMANN  [1976,  Thim.  4.7],  cind  for  the  case  A = an  arbitrary 
commutative  ring,  by  KHARGONEKAR  [1982,  Thm.  3.6].  The  results  we 
have  obtained  are  natural  extensions  of  those  of  Fuhrmann  and 
Khargonekar. 

(4.15)  PROPOSITION.  Let  Q and  Q be  rxr  and  IxZ  invertible 
polynomial  matrices.  A map  ^ : X ^ > X^  _ls  a right  A[z]-module 
isomorphism  if  and  only  if  there  exist  polynomial  [matrices  C,  0,  Y^, 
Yp,  Y^,  and  such  t hat  for  any  x in  Xq 

(4.16)  i|^(x)  = (Cx) 
and  such  that 


(4.17)  CQ  = QD 

(4.18)  CY^  + QY^  = I 


(4.19) 


Y3D  + Y^Q  = I 


□ 
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We  omit  the  proof  of  this  Proposition  since  it  closely  resembles  that 
of  KHARGONEKAR  [1982,  Thm.  3.6].  This  above  proposition  is  used 
later  in  the  dissertation  to  prove  further  results. 

Let  f be  an  input -output  map.  In  general,  there  will  exist 
several  polynomial  matrix  franction  representations  for  W^(z).  We 
now  address  the  following  problem:  If 

W^  = PQ"^R  + S = P Q"^  R + S 

are  two  polynomial  matrix  representations  for  W^  what  conditions 

A A A A 

must  P,  Q,  R,  S,  P,  Q,  R,  and  S satisfy  in  order  that  the  Fuhrmann 
realizations  Z(P,  Q,  R,  s)  and  E(P,  Q,  r\  s)  be  A-isomorphic? 

This  problem  is  closely  related  to  the  problem  of  strict  system 
equivalence  (see  ROSENBROCK  [1970]  and  FUHRMANN  [1977]  for  strict 
system  equivalence  for  systems  over  fields).  Our  results  closely 

resemble  the  results  obtained  by  FUHRMANN  [1977,  Thm.  4.1]  and  by 

KHARGONEKAR  [1982,  Thm.  4.3]  (for  systems  over  rings).  We  state 
without  proof  the 

(4.20)  THEOREM.  Let  f be  an  input/output  map  whose  associated 
transfer  matrix 


W^(z)  = PQ'^R  = p Q'^R 
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Then,  the  Fuhrmann  realizations  Z(P,  Q,  R)  and  l{P,  Q,  R)  are 
A-isomorphic  if  and  only  if  there  exist  polynoinial  matrices  C,  T, 
such  that 


(4.21) 


+ QYg  = I , Y3D  + Y^Q  = I 


□ 


Essentially  equations  (4.21)  oTthe  above  theorem  state  that  the 
system  matrices  (in  the  sense  of  ROSENBROCK  [1970])  must  be 

polynomially  equivalent  for  their  respective  natural  realizations  to 
be  A-isomorphic. 

In  the  next  chapter  we  continue  developing  the  polynomial  model 
theory  and  investigate  Rezout  polynomial  factorizations  of  transfer- 
function  matrices. 


CHAPTER  FIVE 

POLYNOMIAL  FACTORIZATIONS  OF 
TRANSFER-FUNCTION  MATRICES 

In  this  chapter,  we  shall  consider  the  issues  of  existence  and 
computation  of  polynomial  matrix  fraction  representations . Chapter 
Six  explores  application  of  the  theory  developed  here  to  the  design 
of  feedback  control  systems.  The  controllers  designed  will  be 
specified  in  terms  of  a polynomial  matrix  fraction  rep-resentation  of 
the  transfer-function  matrix  and  can  be  implemented  using  the 
polynomial  realization  theory  described  in  the  previous  chapter.  We 
would  like  to  emphasize  that  all  of  our  results  are  constructive. 

Many  frequency -domain  methods  used  to  design  controllers  for 
time-invariant  systems  begin  with  polynomial  factorizations  of  the 
plant  transfer-function  matrix  G(s)  of  the  form 

G(s)  = Q-l(s)R(s),  with 

(5.1) 

Q(s)Yi(s)  + R(s)Y2(s)  = I . 

See  for  example  the  books  of  ROSENBROCK  [1970]  and  WOLOVICH  [1974]. 
Here,  q,  R,  Y^,  and  Y^  are  polynomial  matrices,  and  (5.1)  is  referred 
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3 l6ft-Bezout  polyriorni a1  factorization.  For  time-invariant 
systems,  the  existence  of  these  factorizations  is  not  an  issue;  any 
plant  transfer  function  matrix  always  admits  left-  or  right-Bezout 
polynomial  factorizations.  This  however,  is  not  the  case  with  time- 
varying  systems  (see  Example  (5.9)).  We  therefore  need  to  understand 
what  class  of  time  varying  systems  admit  such  factorizations  before 
we  use  them  to  design  control  systems.  Theorem  (5.2)  answers  just 
this  issue  and  essentially  states  that  a transfer-function  matrix 
admits  Bezout  factorizations  if  and  only  if  the  associated  input/ 
output  map  admits  a canonical  realization. 

(5.2)  THEOREM.  Let  f be  an  input/output  map,  and  let  W^(z)  be  the 
transfer-function  matrix  associated  with  f.  Then,  the  following 
statements  are  equivalent. 

(a)  There  exist  polynomial  matrices  Q,  R,  , and  Yp  over  A[z1 
with  Q invertible  and  such  that 

Wf(z)  = Q"^R  , QYj  + RY^  = I 

— (z ) admits  a right-Bezout  polynomial 

factorization. 

(t*)  There  exist  polynomial  matrices  ?,  Q,  Y-,,  and  Y/|  with 
invertible  and  such  that 


Wf(z)  = PQ"^  , Y3P  + Y^Q^  = I , 


49 


• * llffizJ admits  a 1eft-8ezout  polynomial 

factorization . 

(c)  f admits  a canonicdl  realization. 

PROOr . Tne  proof  of  this  theorem  is  extremely  long  and  rather 

intricate,  and  may  be  found  in  Appendix  B.  We  would  like  to 
emphasize  that  this  proof  is  constructive. 

We  now  present  a systematic  procedure  for  obtaining  the 
polynomial  factorizations  (i.e.,  the  matrices  Q,  R,  Y^,  and  Y2)  of 
part  (a)  of  this  theorem.  An  almost  identical  technique  can  be  used 
to  obtain  the  polynomial  factorizations  of  part  (b)  of  the  above 
theorem. 

Let  Z = (F,  G,  H)  be  (given)  a canonical  realization  of  the 
input/output  map  f. 


_STEP_J_.  Since  the  given  system  Z is  canonical,  it  is  reachable 
in  N steps  for  some  integer  N.  Let  Rf^j  :=  [G  F(aG)  ...  F(oF) 
...  (o' ■-F)(a'^~-^G)]  be  the  N-step  reachability  matrix  for  the 
pair  (F,  G).  Compute  a right-inverse  U for  over  A.  Compute 
the  polynomial  matrices 


0 I zl 

n 0 I 


U 


\ = -(c\) 


0 
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X2  = [I  zl  ...  z^'h]  U . 

Here,  each  of  the  identity  matrices  in  an  nxn  block.  It  will 
happen  that,  with  and  X2  defined  as  above, 

(zl  - F)X^  + GX2  = I 


STEP_H.  Find  matrices  X and  Y over  A such  that 


is  invertible  over  A.  Different  choices  for  the  matrices  X and 
Y will  result  in  unimodularly  related  Bezout  polynomial 
factorizations  (see  Appendix  B). 

STEP  III.  Compute  Nq  :=  Mq^  and  partition  it  as 


where  A is  an  nxn  matrix.  It  will  happen  because  of  the 

observaoi  1 ity  of  the  pair  (F,  H)  (see  Proposition  (8.1))  that 
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the  pair  (A',  B')  is  reachable  in  K steps  for  some  integer  K. 
Determine  (using  Proposition  (2.6))  an  nxp  matrix  S over  A such 

A 

that  A :=  A'  + B'S'  is  a-nilpotent,  i.e.,  for  some  integer  L, 
A(aA)(a^A)  ...  (a‘'~^A)  = 0 . 


STEP  IV.  Define  the  matrix 


W : = 


I 

0 


(a'^S) 

0 


and  compute  the  polynomial  matrix 


V = (2  (-N(,Wz)')N(,  . 


i=0 


Partition  V as 


(5.3)  V = 


'^0 


D -0 


where  Q is  a pxp  matrix. 


^TEP_V.  Define  the  polynomial  matrices 


R :=  DG  , Yj  :=  HXi(zS  + X)  - Y , 
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Y2  = X2(zS  + X) 

Then,  R,  Y^,  Y2  and  Q (as  defined  in  5.3)  determine  the  desired 

A 

polynomial  factorization  of  W^(z)  as  in  Theorem  (5.2). 

We  now  illustrate  this  procedure  using  two  examples. 

(5.4)  EXAMPLE.  Consider  the  linear  time-varying  system  X = (F,  g, 
h)  where  F(k)  = 1 for  all  k in  Z,  and  g = h = a,  where  a(k)  = 1 if  k 
is  even  and  a(k)  = 0 otherwise  (see  Example  (2.5)).  We 
systematically  employ  the  procedure  described  earlier  to  obtain  a 
right-Bezout  polynomial  factorization  for  W^(z)  = h(zl  - F)"^g. 

STEP — The  pair  (F,  g)  is  easily  seen  to  be  reachable  in  two 
steps,  and  R2  = [g  F(ag]  = [a  1-a].  A right-inverse  for  R2  is 
simply  U = [1  1].  Then, 


'0 

f 

= 

-(aR2) 

U 

= a - 1 

_0 

0_ 

X2  = 

[1 

z]  U = 

z + 1 

STEP  II. 

With  X = 

0, 

and  Y 

= 1. 

the  matri 

-F 

X 

-1  0~ 

% 

: = 

= 

Y 
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is  invertible  over  A. 

STEP  III. 

■-1  0~ 

a 1 

We  need  to  find  an  S in  A such  that  (-1  + aS)  is  a-nil potent. 
By  inspection,  it  is  clear  that  S = 1 will  suffice. 


:=  M 


-1 


STEP  IV.  Defining 


W : = 


1 1 
0 0 


we  compute  the  polynomial  matrix  V: 


V = ( E (-N  Wz)')  N 
i =0  ^ ^ 


-1-az 

az^+z 

a+az 

2 

-az  -az+1 

STEP  V.  The  desired  factorization  is  then 
W^(z)  = Q‘^R  , QY^  + RY^  = 1 , 
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where 


Q = az^  + az  - 1 , R = a 


2 

+ z 


□ 


(5.5)  EXAMPLE.  Let  us  consider  an  armature-control  1 ed  dc  motor 
described  by  the  input/output  differential  equation 


Here,  the  input  u(t)  is  the  applied  armature  voltage,  the  output 
0(t)  is  the  angular  position  of  the  motor  shaft,  and  the  time- 
varying  coefficients  a(t)  and  3(t)  are  given  by 

a(t)  = +0.74  + 0.3k(t)  , 3(t)  = 5.56k(t) 

where  k(t)  is  the  (normalized)  effective  motor  torque  "constant". 
The  nominal  value  of  k(t)  is  1,  however  during  operation  k(t)  may 
vary  significantly  due  to  motor  shaft  loading  and  heating  effects. 
We  assume  that  the  time-variation  of  k(t)  is  known  a priori;  for 
example  the  motor  may  be  part  of  a machining  operation  and  k(t)  could 
be  identified  off-line,  say  by  measurements  taken  during  test  runs. 


(5.6) 


+ a(t)  = e(t)u(t) 
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By  sampling  the  continuous-time  system  (5.6)  with  a suitably 
small  sampling  period  T we  obtain  the  following  time-varying  sampled- 
data  state  model  Z =(F,  G,  H)  that  approximates  the  motor 
characteristics: 


■ e(kT+T)“ 

■ 1 T ~ 

■0(kT)  ■ 

" 6T^/2" 

= 

+ 

_o)(kT+T)_ 

_ 0 1-«T_ 

_o)(kT)  _ 

- ST 

u(kT) 


> 


(5.7) 


6(kT)  = [l 


•0(kT) 

“(kT), 


Here  oj(kT)  is  angular  velocity  fo  the  motor  shaft.  It  can  be  easily 
verified  that  Z = (F,  G,  H)  is  canonical.  We  wish  to  obtain  a right- 

Bezout  polynomial  factorization  for  the  transfer-function  matrix 

~ -1 

Wj,(z)  = H(zl  - F)“  G.  We  shall  later  use  this  factorization  (see 
Example  (6.7))  to  design  a dead-beat  controller  for  the  motor.  We 
systematically  employ  the  procedure  described  earlier. 


STEP  I.  It  can  easily  be  shown  that 
(zl  - F)Xj  + GX2  = I 


where 
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-T/2 

2 

T74 

-1 

T/2 

} C-z  + (1  - ctT) 


X 


STEP  II.  With  X'  = [0 


1/T],  Y = 0 


-F 

X 

~-l 

-T 

= 

0 

-1+aT 

H 

Y 

1 

0 

is  invertible  over  A. 


STEP  III. 


N,^  :=  M, 


-1 


0 

-1/T 

-1+aT 


0 

0 

T 


I 

X 


, the  matrix 


0 

1/T 

0 


1 

1/T 
1+aT  . 


Notice  that  the  matrix  A :=  the  upper  2x2  block  of  Nq  is 
already  o-nil potent.  Hence  we  may  choose  S = [0  0]'. 
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STEP  IV.  Defining 


I 

(o-'s)' 

1 

0 

0 

w = 

= 

0 

1 

0 

_0 

0 

_0 

0 

0_ 

we  compute 


V : 


N 


0 


0 

-1/T 


z+aT-1 


0 1 

0 1/T(z-1) 

T -z^+(2-aT)z+aT-l 


STEP  V.  The  desired  factorization  is  then 


Wj.(z)  = Q-^R 


QYj  + RY2  = 1 , 


where 


(5.8) 


2 2 

Q = Z^  + (aT-2)z  + (1  - aT)  , R = + -^(aT+1)  , 

^ ^ ’ ^2  = XT  • 


aT‘ 


Recall  that  X 


2 


+ T. 


□ 
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We  would  like  to  remark  that  not  all  time -varying  systems  admit 
canonical  realizations.  This  is  demonstrated  by  the  following 
example. 

(5.9)  EXAMPLE.  Consider  the  linear  time  varying  system  L = 
(F,  g,  h)  where  F(k)  = h(k)  = 1 for  all  k in  Z,  and  g = A = the  unit 
pulse  concentrated  at  the  origin.  Let  fj.  be  the  input/output  map 
associated  with  E.  We  will  show  that  fj.  does  not  admit  a canonical 

A 

realization,  and  therefore,  by  Theorem  (5.2)  W^(Z)  does  not  admit 

left-  or  right-  polynomial  factorizations. 

Suppose  fj.  did  admit  a canonical  realization.  Then,  by  Theorem 
(5.2),  there  exist  polynomial  matrices  P and  Q with  Q invertible  -and 
such  that 


(5.10)  h(zl  - F)-lg  = (z  - I)'!  A = PQ'l  . 
n 

Let  0 = E z^Q. . From  (5.10)  it  follows  that  (z  - 1)~^A  Q = P is  a 
i=0  ^ 

polynomial.  The  coefficient  of  z"^  above  must  be  zero: 

0 = AQg  + (aA)Q^  + ...  + (o%)Q^  . 

Multiplying  this  equation  on  the  left  by  (a^A)  we  conclude  that 
(a^A)Q^.  = 0 for  i = 0,  1,  ...,  n.  Thus  AQ(z)  = 0 which  is 

impossible  since  Q is  assumed  to  be  invertible.  CH 
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The  central  theorem  of  this  chapter,  Theorem  (5.2),  is  stated 
for  canonical  linear  time-varying  systems.  Under  the  weaker 
hypothesis  of  reachability  we  can  obtain  the  following  theorem  (which 
is  actually  an  intermediate  result  used  in  proving  Theorem  (5.2)), 
the  proof  of  which  may  be  found  in  Appendix  B. 

(5.11)  THEOREM.  Let  T.  = (F,  G,  H)  be  a linear  time-varying  system 
that  is  reachable  in  N steps.  Then,  there  exist  polynomial  matrices 
P and  Q with  Q invertible  such  that 


Wj.(z)  = H(zl  - F)"^G  = PQ~^ 

and  the  Fuhrmann  realization  £(P,  Q,  I)  is  A-isomorphic  to 
T.  = (F,  G,  H) . I — I 


We  would  like  to  remark  in  closing  that  the  systematic  procedure 
given  in  this  section  for  computing  Bezout  polynomial  factorizations 
for  time-varying  systems  carries  through  with  minor  modifications  for 
systems  over  a principal  ideal  domain.  The  interested  reader  may 
refer  to  POOLLA  and  KHARGONEKAR  [1983]  for  more  details. 

In  the  next  chapter,  we  apply  the  theory  developed  here  and  in 


Chapter  4 to  the  design  of  feedback  control  systems. 


CHAPTER  SIX 

APPLICATIONS  OF  THE  POLYNOMIAL  THEORY  TO  FEEDBACK  CONTROL 

We  will  now  explore  applications  of  the  theory  developed  in 
earlier  chapters  to  the  design  of  feedback  control  systems.  In 

particular,  we  will  show  that  via  dynamic  output  feedback  it  is 
possible  to  "coefficient  assign"  canonical  systems.  We  would  like  to 
stress  that  our  design  techniques  are  constructive . The  controllers 
we  obtain  will  be  specified  in  terms  of  a polynomial  matrix  fraction 
representation  of  their  transfer-functions,  and  can  be  then 
implemented  (realized)  using  the  polynomial  realization  theory  of 
Chapter  Four. 

Let  us  consider  a canoni cal  linear  time-varying  system 

A 

Z = (F,  G,  H)  over  A.  From  Theorem  (5.2)  it  follows  that  W„(z)  = 

la 

H(zl  - F)"^G  admits  a right-Bezout  polynomial  factorization: 

(6.1)  Wj.(z)  = Q"^R  , QY^  + RY2  = I . 

Consider  now  the  feedback  control  system  shown  below. 
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V(z) 


E(z) 


■> 


U(z) 


■> 


■> 


Y(z) 


Here  is  the  controller  (to  be  designed).  Suppose  E^  is  specified 
in  terms  of  the  following  collection  of  input/output  difference 
equations : 


U(z)  = P^C(z) 

(6.2) 

0^5(z)  = R^E(z)  . 

This  is  the  so-called  Rosenbrock  representation  (see  discussion  at 

the  end  of  Chapter  Three),  and  here  P^,  0^,  and  R^.  are  polynomial 

matrices  over  A[z]  and  is  assumed  to  be  invertible.  Then  the 

controller  transfer-function  matrix  has  the  form  (z)  = P 0“  R . 

E^ ' ' c c c 
c 

We  can  then  represent  the  closed-loop  system  shown  above  by  the 
equati ons 


(6.3) 


' Q 

-Rp; 

"y(z)“ 

" 0 ' 

-R 
L c 

- 

J(z)_ 

R 

L c J 
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Thus,  we  see  that  the  closed-loop  system  dynamics  are  determined  by 
the  inverse  of  the  matrix 


tt(z)  = 


In  particular,  the  closed  loop  system  is  internally  uniformly 
asymptotically  stable  if  and  only  if  is  a stable  Laurent 

series  (see  Chapter  Seven). 

We  now  address  the  following  question:  To  what  extent  can  the 
matrix  ir'^(z)  be  "assigned"  by  selecting  P^,,  Q^,  and  R^,  i .e.  by 
designing  the  controller.  Recall  that  we  must  constrain  to  be 
invertible  and  W^.  (z)  = to  be  proper  in  order  to  realize 

the  controller.  In  the  time-invariant  case  this  problem  has  been 

extensively  studied.  It  has  been  shown  that  (see  for  example  EMRE 

and  KHARGONEKAR  [1982])  with  R^.  = I,  it  is  possible  to  arbitrarily 
assign  the  coefficients  of  the  polynomial 

detir(z)  = det(QQ^  - RP^) 

For  the  time-invariant  case,  since  the  roots  of  the  equation 

detir(z)  = 0 are  the  closed-loop  system  poles,  this  assignability 
result  implies  that  it  is  possible  to  arbitrarily  alter  the  closed- 
loop  system  dynamics  by  appropriately  designing  the  controller  (i.e. 
selecting  P^  and  Q^.) . 
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In  the  time-varying  case,  the  matrix  ir(z)  is  defined  over  the 
noncommutative  ring  A[z]  and  there  is  no  known  definition  for  the 
determinant  of  a matrix  over  a non-commutati ve  ring.  Thus,  we  cannot 
speak  of  assigning  the  coefficients  of  detir(z).  However,  we  can 
still  consider  the  problem  posed  earlier;  to  what  extent  can  tt~^(z) 
be  assigned  by  choosing  W^,  (z).  In  this  regard,  we  have  the 
following  result: 


(6.4)  THEOREM.  Let  Z = (F,  G,  H)  be  a canonical  linear  time-varying 
m-input  p-output  system  over  A,  and  let 


Wj,(z)  = H(zl  - F)'^G  = Q'^R  , QY^  + RY^  = I 


be  any  pxp  monic  polynomial 

matrix 

wi  th 

deg(<t>)  > 2 deg(Q)  + 

max{deg( Y^) ,deg( Y2)-1} . Then, 

there  exists 

a controller  E with 

proper  transfer  function  matrix 

W.  (c) 

= P Q" 

^R„  where  P^,  0^  and 

“c 

R^  are  polynomial  matrices,  and  there 

exi  St 

unimodular  polynomial 

matrices  T^(z)  and  Tg(z)  such  that 


= Ti(z) 


I 

0 


T2(z) 


64 


PROOF.  Since  Q is  an  invertible  polynomial  matrix,  it  follows  from 
Proposition  (3.2)  that  there  exists  a polynomial  matrix  T such  that 
M :=  TQ  is  monic  and  deg(M)  = deg(Q).  We  now  divide  (}>  on  the  left  by 
TQ  to  obtain 

(6.5)  (|)  = MjTQ  + N 

where  deg(N)  < deg(Q).  Define  the  polynomial  matrices  P^,  Q^,  and 
by 

(6-6)  Pc  :=  Y2  , Qc  :=  M^T  + NY^  , R^  :=  -N  . 

We  shall  show  that  the  desired  controller  has  transfer  function 
matrix  “ ^c^c  *^c  '^bere  P^,  Q^,  and  R^  are  defined  by 

(6.6) .  Before  we  do  this,  we  must  first  ensure  that  (i ) is 

invertible  and  that  (ii)  PcQc^Rc  is  proper. 

To  prove  (i),  notice  that 

deg(T-^)  = deg(QM“l)  = deg(Q)  - deg(M)  = 0 , 

deg(Mj)  = deg((t>)  - deg(M) 

> deg(Q)  + deg(Y^) 

It  then  follows  that 

Qc  = MiT(I  + T-^M];1nYj)  , 


and  that 
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deg(T-lMilNYi)  < deg(T-l)  - deg(M^)  + deg(N)  + deg(Yi) 

< 0 - deg(Q)  - deg(Y^)  + deg(N)  + deg(Y^) 

< deg(N)  - deg(Q) 

< 0 . 

Thus,  I + is  monic,  and  hence  is  invertible. 

We  now  prove  (ii).  Notice  that  we  can  write 


Pc^c^^c  " '^2^^  ^ T“^M‘^NY^)‘^T'^M^^N  . 


However,  since  deg(T-l)  = deg(I  + T-1m];1nYi)-1  = 0,  it  follows  that 
deg(P^Q-lR^)  < deg(Y2)  + deg(N)  - deg(M^) 

< 0 . 

Thus,  the  controller  transfer-function  matrix  W_  (z)  = P 0"^R  is 

' ' c c c 

well  defined. 

With  the  definitions  (6.6)  for  P^.,  Q^,  and  R^,  it  can  now  be 
mechanically  verified  that 


Q 


I 


I 0 
4> 


=:  T^(z) 


0 


T2(z)T3(z) 
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The  above  equation  together  with  the  observation  that  T2, 
and  T3  are  unimodular  polynomial  matrices  completes  the  proof. 

Our  technique  in  the  above  proof  closely  resembles  that  of 
KHARGONEKAR  and  OZGULER  [in  press]. 

(6.6)  REMARK.  We  could,  for  instance,  choose  (|)(z)  = z''l  with  r a 
sufficiently  large  integer  to  satisfy  the  degree  constraint  on  ({>.  In 
this  event,  the  controller  designed  as  in  the  above  theorem  will 
result  in  a "dead-beat"  closed-loop  system  response.  In  other  words, 
with  this  controller,  the  closed-loop  system  response  resulting  from 
any  initial  state  will  become  zero  after  a finite  number  of  steps 
(assuming  the  external  input  V(z)  = 0).  □ 

We  now  present  a systematic  procedure  for  designing  a controller 
= (F^,  G(.,  H^)  to  yield  a desired  closed-loop  system  response 
(specified  by  the  polynomial  matrix  (t>).  This  procedure  is 
essentially  distilled  from  the  proof  of  the  previous  theorem  and  is 
exhibited  explicitly  for  the  reader's  convenience. 

STEP  I . Given  a canonical  linear  time-varying  system  E = (F,  G,  H), 
first  obtain  a right-Bezout  polynomial  factorization  for  the 
transfer-function  matrix  as 

Wj.(z)  = Q'^R  , QY^  + RY2  = I 
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This  can  be  done  by  systematically  employing  the  procedure  described 
in  the  previous  chapter. 

STEP  II.  Let  (J)(z)  be  a (given)  monic  pxp  polynomial  matrix  which 
represents  the  desired  closed-loop  system  response  as  in  the 
statement  of  Theorem  (6.4).  Recall  that  there  is  a contraint  on  (j)  : 
deg((j))  > 2 deg(Q)  + max{deg(Y^),  deg(Y2)  - 1).  Left -divide  <t>  by  0 to 
obtai n 


(|)  = XQ  + N 
with  deg(N)  < deg(Q) . 

STEP  III.  Compute  the  polynomial  matrices 

Pc  :=  Y2  , Qc  :=  (X  + NY^)  , R^  :=  -N  . 

A 

Then,  the  controller  transfer-function  matrix  is  given  by  W^.  (z)  = 

STEP  IV.  Realize  the  controller  by  a state  model  = (F^.,G^,H^.) . 
This  can  be  done  using  the  polynomial  realization  theory  of 
Chapter  4.  If  Q is  monic,  it  will  happen  that  Q^,  is  also  monic,  in 
which  case  explicit  formulae  for  the  state  model  = (F^,G^,H^,J^.) 
are  given  by  (4.11) . 
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We  conclude  this  chapter  with  an  example  illustrating  the  design 
procedure  outlined  above. 

(6.7)  EXAMPLE.  Consider  again  the  armature-control  led  dc  motor 
described  in  Example  (5.5),  with  a time-varying  motor  torque 
"constant".  The  time-variation  could  be  due  to  motor-shaft  loading 
or  heating  effects.  We  shall  use  the  systematic  procedure  described 
above  to  design  a "dead-beat"  controller  for  the  motor. 

STEP  I.  The  time-varying  sampled-data  state  model  Z = (F,  G,  H)  of 
the  motor  is  specified  by  Equations  (5.7).  We  have  already,  in 
Example  (5.5)  computed  a right-Bezout  polynomial  factorization  for 
Wj.(z) , and  obtained 


STEP  II.  Since  we  wish  to  design  a "dead-beat"  controller  for  E we 


W^(z)  = Q’^R  , QY^  + RY^ 


where 


Q = z^  + (aT  - 2)z  + (1  - aT) 


^1  ■ 4X  ’ ^2  " 20^  20^ XT 


choose  (j)(z)  = z^  [see  Remark  (6.6)].  Notice  that  <J)  satisfies  the 


69 


degree  constraint  deg(<|.)  > 2deg(Q)  + max{degYj,  degY2  - 1}.  Left- 
dividing  <t>  by  Q as  (}i  = XQ  + N we  obtain 

X(z)  = z2  - z(a-l(aT  - 2))  + c , 

N(z)  = [a-l(aT  - 2)  a'^d  - aT)  - (aT  - 2)c]z  - (1  - ctT)c  , 

c = a"2(aT  - 2)a"l(aT  - 2)  - a“l(l  - aT) 

STEP  III.  The  controller  transfer-function  matrix  is  then  given  by 

Wj.  (z)  = - Y^CX  + NY^)"^N 

and  can  readily  be  realized  using  the  formulae  (4.11)  since  in  this 
case,  Q(z)  is  monic.  We  leave  the  details  to  the  interested  reader.lZ] 

In  this  first  part  of  the  dissertation  we  have  examined  in 
3 polynomial  theory  for  linear  time-varying  systems.  For  the 
remainder  of  this  dissertation  we  turn  our  attention  to  the  study  of 
linear  time-varying  systems  based  on  the  ring  of  stabl e-proper 
rational  functions  in  Jl“(Z)[[z"^]]. 


CHAPTER  SEVEN 
STABILITY 

In  this  chapter,  we  briefly  review  some  important  well-known 
concepts  dealing  with  the  stability  of  linear  time-varying  systems. 
We  also  relate  some  of  these  concepts  to  the  transfer-function  theory 
based  on  skew  rings  developed  in  this  dissertation.  Some  of  the 
material  in  this  chapter  is  condensed  from  Section  4 of  KAMEN, 
KHARGONEKAR  and  POOLLA  [1984].  The  reader  is  referred  to  this  paper 
for  further  details. 

In  the  past,  a good  deal  of  research  on  the  stability  of  both 
continuous-  and  discrete-time  linear  time-varying  systems  has  been 
done.  For  example,  the  interested  reader  may  consult  articles  by 

WILLEMS  [1970],  ANDERSON  and  MOORE  [1969,  1981a]  (for  Liapunov 

stability  theory),  CESARI  [1963],  STARZINSKII  [1955]  (for  stability 
of  special  classes  of  time-varying  systems  such  as  periodic  systems), 

FREEDMAN  and  ZAMES  [1968]  (for  stability  of  si owly-varyi ng  systems), 
etc. 

For  any  vector  x in  R"  let  iixil  denote  the  Euclidean  norm  of  x. 

For  an  mxn  matrix  M over  R,  let  HMH  :=  sup  dMxII/HxH.  For  an  nxm 

matrix  N over  A let  the  norm  of  N be  defined  by  IIN»  :=  sup  HN(t)il. 

t 
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Let  f : a'J  -»■  aP  be  an  input/output  map.  We  shall  say  that  f is 

bounded  input-bounded  output  (BIBO)  stable  if  and  only  if  for  any 

bounded  input  sequence,  i.e.,  for  any  u over  ^(Z)^.  the 

corresponding  output  sequence  y = f(u)  is  bounded,  i.e.,  y is  over 
00 

A (Z)^.  Let  I.  = (F,  G,  H,  J)  be  any  (fixed)  realization  of  f. 
Consider  the  free  behaviour  of  the  system  Z described  by  the  vector 
difference  equation 

(7.1)  x(k  + 1)  = F(k)x(k)  . 

The  system  I is  said  to  be  internally  uniformly  asymptotically 

^ (u.a.s.)  if  and  only  if  for  every  real  number  e > 0,  there 

exists  a positive  integer  such  that  for  any  initial  time  t^  in  Z 
and  any  initial  state  x(tg)  with  nx(t^)ii  < 1,  we  have  that 

iix(tQ  + i)n  < e for  all  i > N^.  Here,  x(tg  + i)  is  the  solution 
of  (7.1)  at  time  t^  + i starting  from  initial  state  x(t^). 

Let  P(z)  be  an  nxm  matrix  over  the  skew  ring  A((z“^)),  i.e.  P(z) 
is  a matrix  Laurent  series  of  the  form 

oo 

P(z)  = E z"‘'p. 
i=-N  ^ 

Then,  P(z)  is  said  to  be  stable  if  and  only  if 
n P^.  a ->■  0 as  i -»■  <» 
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In  terms  of  this  notion,  we  can  characterize  internal  stability  as 
follows  (see  GREEN  and  KAMEN  [1983]  and  KAMEN,  KHARGONEKAR  and  POOLLA 
[1984,  Prop.  (4.4)]). 

(7.2)  PROPOSITION.  Let  I = (F,  G,  H)  be  a linear  time -varying 
system  over  A.  Then,  Z is  internally  u.a.s.  if  and  only  if  (zI-F)~^ 
is  a stable  matrix  power  series.  □ 

Note  the  similarity  this  result  bears  to  the  time-invariant 
theory.  Let  f be  an  input/output  map,  and  let  Z = (F,  G,  H,  J)  over 
A be  any  (fixed)  realization  of  f.  In  contrast  with  the  time- 
invariant  case,  internal  stability  of  Z does  no^  in  general  imply 
BIBO  stability  of  f.  However,  for  the  class  of  bounded  linear  time- 
varying  systems,  we  have  the  following  result  (the  proof  is  omitted 
on  account  of  its  relative  ease). 

(7.3)  PROPOSITION.  Let  f be  an  input/output  map  and  let  Z = 
lf.»  0,  H,  J)  over  i (Z)  be  any  (fixed)  realization  of  f.  Suppose  Z 
is  internally  u.a.s.  Then  f is  BIBO  stable.  Q] 

The  following  examples  illustrate  phenomena  peculiar  to  time- 
varying  systems  with  unbounded  coefficients. 


(7.4)  EXAMPLE.  Consider  the  time-varying  system  Z = (F,  G,  H)  over 
A where  F = 1/2,  G = 1,  and 
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From  Proposition  (7.2)  it  is  evident  that  E is  internally  u.a.s. 
This  system,  however,  is  not  BIBO  stable  because  the  output  y(k) 
resulting  from  the  input  u = A = unit  impulse  at  the  origin  grows 
without  bound.  Notice  that  Proposition  (7.3)  does  not  apply  here 
because  H is  not  over  t°°(Z). 


(7.5)  EXAMPLE.  Consider  the  linear  time-varying  system  E = 
(F,  0,  H)  over  A defined  by  G = H = 1,  and 


F(k) 


iO,  k<0 

1,  k>0 


Define  a time-function  T over  A by 


T(k) 


k<0 

k>0 


Notice  that  T is  not  an  element  of  I (Z)  and  that  T has  an  inverse 
T"^  over  A.  Consider  the  system  E = (F,  G,  R)  over  A where 


74 


F = (a‘^T‘^)FT  , R = HT 


5 = (a"^r^)G 


The  systems  E and  E are  algebraically  A-isomorphic  (see 
Chapter  2).  However,  E is  not  internally  u.a.s.  while  E 2_s 
internally  u.a.s.  This  is  because  the  transformati on  R is  not 
Liapunov  (recall  that  an  nxn  matrix  T is  a Liapunov  transformati on  if 


Let  f be  an  input-output  map  with  associated  transfer-function 

A 

matrix  W^(z).  Suppose  there  exist  polynomial  matrices  P,  Q,  and  R 

A 

of  appropriate  sizes  with  Q monic  and  such  that  W^(z)  admits  the 
polynomial  representation 


In  Chapter  Four,  we  have  described  in  detail  a technique  to  obtain 
the  (natural)  Fuhrmann  realization  E(P,  Q,  R)  = (F,  G,  H,  J)  of  f 
associated  with  the  particular  polynomial  representation  (7.6)  (the 
explicit  formulae  for  F,  G,  H and  J are  given  in  equation  (4.11)). 
The  following  theorem  characterizes  internal  stability  of  the 
Fuhrmann  realization  in  terms  of  the  polynomial  representation  (7.6). 


1 oo 

and  only  if  T and  T“^  are  over  £ (Z)). 


□ 


(7.6) 


W^(z)  = PQ'^R 
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(7.7)  THEOREM.  Let  f be  an  input/output  map  whose  associated 

A 

transfer  function  matrix  W^(z)  admits  the  polynomial  representation 
W^(z)  = PQ'^R.  Further  assume  that  P,  Q,  and  R are  polynomial 
matrices  over  Z (Z)  [z],  with  Q monic.  Then,  the  Fuhrmann 
realization  E (P,  Q,  R)  = (F,  G,  H,  J)  is  internally  u.a.s.  if  and 
only  if  Q~^(z)  is  a stable  matrix  Laurent  series. 


PROOF.  Since  the  Furhmann  realization  associated  with  the  polynomial 
factorization  W^(z)  = H(zI-F)"^G  + J is  itself  2 = (F,  G,  H,  J),  it 
follows  that  in  particular  Xg  and  X^j_p  must  be  isomorphic  as  right 

oo 

IL  (Z)[z]-modules.  This,  by  Proposition  (4.15)  implies  that  there 
exist  polynomial  matrices  C,  D,  and  ^2  ^”(Z)[z]  such  that 

CQ  = (zI-F)D  , + YgO  = I 

From  these  equations,  we  can  write 

Q'^  = YiDQ"^  + Y^ 

= Y^(zl  - F)"^C  + Y2  . 

It  is  now  clear  from  the  last  equality  that  if  (zl  - F)“^  is  stable 
then  0“^  is  stable.  The  converse  follows  from  an  identical  argument. 

□ 
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In  subsequent  chapters,  we  shall  deal  only  with  linear  time- 
varying  systems  with  bounded  coefficients,  i.e.,  defined  over  the 
difference  subring  Z (Z)CA.  We  shall  also  require  that  any 
controller  we  design  be  over  i (Z).  These  are  physically  reasonable 
constraints  since  most  time-varying  plants  that  arise  in  practice 
have  bounded  time  variation,  and,  implementation  of  controllers  with 
unbounded  coefficients  would  be  numerically  ill-conditioned. 
(Technically,  these  constraints  substantially  complicate  proofs  and 
make  results  harder  to  obtain). 


/ 


CHAPTER  EIGHT 

STABILITY  AND  ASYCONTROLLABILITY 

In  this  chapter,  we  introduce  the  key  notion  of 
asycontrol 1 abi 1 ity , which  is  closely  related  to  being  able  to 

stabilize  a linear  time-varying  system  by  dynami c state  feedback. 
ANDERSON  and  MOORE  [1981a]  have  defined  the  notion  of  stabi 1 i zabi 1 ity 
which  is  equivalent  to  being  able  to  stabilize  a linear  time-varying 
system  by  nondynami c (i.e.,  memoryless)  state  feedback.  One  of  our 
main  results  in  this  chapter.  Theorem  (8.9),  shows  the  equivalence  of 
stabi li zabi 1 ity  and  asycontrol lability.  striking  conclusion  of 

this  theorem  is  that  dynami cs  in  state  feedback  buy  nothing  extra  as 
far  as  the  problem  of  stabilization  is  concerned. 

Let  us  consider  a linear  time-varying  system  l = (F,  G,  H,  J) 

defined  over  the  difference  subring  Ji”(Z)  A.  Recall  that  by 

Corollary  (4.14),  Z is  reachable  in  N steps  if  and  only  if  there 

exist  polynomial  matrices  Y^  and  Y2  such  that 

(8.1)  (zl  - F)Yi  + GY2  = I 

Motivated  by  equation  (8.1)  we  have  the  following  key 
definition: 
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(8.2)  DEFINITION.  The  system  Z is  said  to  be  asycontroHable  if  and 
only  if  there  exist  stable  matrix  Laurent  series  and  Yq  such  that 

(zl  - F)Yi  + GY2  = I . □ 

This  technical  definition  has  a rather  complex  precise  system- 
theoretic  interpretation  (see  Appendix  C and  also  KHARGONEKAR  and 
POOLLA  [1984b])  in  terms  of  stabilizing  I via  an  open-loop  control 
law.  Roughly  speaking,  a system  is  asycontrol 1 able  if  and  only  if  it 
can  be  driven  to  zero  "final"  state  asymptotically,  using  uni formly 
bounded  input  sequences  along  uni formly  bounded  state  trajectori es . 
In  particular,  this  implies  that  if  E can  be  stabilized  using  a 
dynamic  state-feedback  controller,  then  E is  ascontrollable.  The 
phrase  asycontrol 1 abi  1 i ty  (from  asymptotically  control lable)  is 
borrowed  from  KHARGONEKAR  and  SONTAG  [1982]  for  time-invariant 
systems  over  rings. 

We  shall  also  need  the  following  definition.  A system  E = 
(F,  G,  H,  J)  over  l (Z)  is  called  ri ght-rational  asycontrol 1 able  if 
and  only  if  there  exist  polynomi al  matrices  A,  B,  and  S over 

00  4 

I (Z)[z]  with  S invertible  and  a stable  Laurent  series  and  such 
that 


(zl  - F)AS"1  + GBS-1 


I 
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We  now  relate  right-rational  asycontrol lability  to  the  existence  of  a 
stabilizing  controller  with  the  following  result. 

(8.3)  THEOREM.  Let  I = (F,  G,  I)  over  £ (Z)  be  a linear  time- 
varying  system  that  is  right-rational  asycontrollable.  Then,  there 
exists  a controller  = (F^,  H^)  over  &°°(Z)  such  that  the 

closed  loop  system  defined  by 


X 

Pf  GH  1 

X 

r-  “ 

G 

c 

(k  + 1)  = 

(k) 

(k)  + 

X 

G.  F 

L c c J 

0 _ 

is  internally  uniformly  asymptotically  stable  (u.a.s.). 

SKETCH  OF  PROOF.  We  would  like  to  remark  that  the  proof  of  this 
result  is  constructive  enabling  one  to  explicitly  compute  the 
controller  using  operations  in  the  skew-ring  A[z].  This  proof 
closely  resembles  the  proof  of  Theorem  (6.4).  We  outline  the  basic 
elements  of  the  construction. 

(i)  Define  K :=  max{deg(AS“^) , deg(BS"^)  - 1).  Pick  any  moni c 
polynomial  matrix  over  £*(Z)[z]  with  deg(Sj^)  > K + 2 
and  such  that  is  a stable  matrix  Laurent  series. 


(ii) 


Divide  on  the  left  by  (zl  - F)  to  obtain 
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S-^  = M(zl  - F)  + N 
with  deg(N)  = 0. 

(iii)  Define  the  polynomial  matrices 

P(.  :=  B , :=  MS  + NA  , :=  -N  . 

Then,  the  controller  transfer  function  is  function  is 

given  by  W^.  (z)  = • 

c 

(iv)  Implement  the  controller  using  the  polynomial  realization 
theory  of  Chapter  4.  Explicit  formulae  for  the  controller 
in  terms  of  the  coefficient  matrices  of  P^.,  Q^,  and  are 
given  (for  the  case  when  Q^.  is  monic)  by  equations  (4.11). 

□ 

Let  2 = (F,  G,  H)  be  a linear  time-varying  system  over  Jt"(Z). 
ANDERSON  and  MOORE  [1981a]  have  defined  a notion  of  stabilizability 
which  intuitively  corresponds  to  requiring  that  unstable  modes  be 
controllable.  The  authors  then  show  that  this  notion  is  equivalent 
to  the  existence  of  a stabilizing  memory less  state  feedback  law 
u(k)  = -L(k)  x(k).  We  shall  (perversely)  take  this  to  be  the 

definition  of  stabilizability.  In  terms  of  our  skew-ring  framework, 
we  phrase  this  as  follows: 
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(8.4)  DEFINITION.  Let  T.  = (F,  G,  H)  over  jt”(Z)  be  a linear  time- 

varying  system.  Then,  Z is  said  to  be  stabilizable  if  and  only  if 

there  exists  an  mxn  matrix  L over  such  that  (zl  - F + GL)~^ 

is  a stable  matrix  power  series. 

□ 

One  can  similarly  define  the  dual  notion  of  detectabi 1 i ty . 

Let  Z = (F,  G,  H)  over  a”(Z)  be  a stabilizable  linear  time- 
varying  system.  Let  L be  an  mxn  feedback  matrix  as  in  Definition 

(8.4) .  Notice  that  we  can  write 

(zl  - F)Y]^  + GY2  = I 

where  Y^  = (zl  - F + GL)~^  ^2  ~ *-^l*  Since  Y^  and  Y2  are  stable 

Laurent  series,  it  follows  from  Definition  (8.1)  that  Z is 
asycontrol lable.  Thus,  stabilizab^ility  impl i es  asycontrol  1 abi  1 i ty . 
The  much  more  difficult  converse  is  also  true. 

(8.5)  THEOREM.  Let  Z = (F,  G,  H,  J)  be  a linear  time-varying  system 
over  I (Z).  Then,  Z is  asycontrol  1 abl e if  and  only  if  Z is 
stabilizable. 

PROOF.  The  proof  of  this  theorem  is  extremely  long  and  may  be  found 
in  Appendix  C.  We  would  like  to  remark  that  the  essential  technical 
difficulty  in  the  proof  is  ensuring  that  the  feedback  matrix  L is 
over  t”(Z).  n 
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We  also  immediately  obtain  the  following: 

(8.6)  COROLLARY.  A system  H is  right-rational  asycontrollable  if 
and  only  if  it  is  asycontrollable. 

PROOF.  We  have  the  following  sequence  of  implications  proving  our 
cl  aim. 


stabil izabil ity 


(a) 

— ^ right-rational  asycontrol labil ity 

(b) 

V 

asycontrol labil  ity 


Here,  (a)  follows  from  the  discussion  preceding  Theorem  (8.5),  (b) 
follows  trivially  from  the  definition  of  right-rational 
asycontrollability,  and  (c)  is  a restatement  of  Theorem  (8.5).  CH 

Recall  (from  the  discussion  following  Definition  (8.2))  that  if 
a linear  time-varying  system  Z can  be  stabilized  by  a dynamic  state- 
feedback  controller,  then  Z is  asycontrollable.  This  observation 
together  with  the  Definition  (8.4)  of  stabilizability  immediately 
offers  the  following  surprising  conclusion. 


83 


(8.7)  COROLLARY.  If  a linear  time-varying  system  I can  be 

stabilized  using  dynamic  state-feedback,  then  Z can  also  be 
stabilized  using  memoryless  state-feedback.  [ | 

(8.8)  REMARK.  Indeed  the  above  corollary  is  not  (necessarily) 

expected,  because  there  are  classes  of  systems  (for  example  delay 
systems,  see  KAMEN  [1982,  Ex.  3,  p.  371])  for  which  it  is  not  true. 
For  time-invariant  systems  over  a field,  it  is  a well-known  fact  that 
dynamic  state  feedback  is  equivalent  to  memoryless  state  feedback  as 
far  as  the  problem  of  stabilization  is  concerned.  The  proof  of  this 
fact  (see  any  modern  text  on  control  theory)  relies  heavily  on  the 
Kalman  canonical  decomposition.  No  such  decomposition  exists  for 
time-varying  systems  because  the  "dimension"  of  the  reachable  space 
could  depend  on  time.  QJ 

Since  stabi 1 izabi 1 ity  and  asycontrol lability  are  equivalent 
notions,  we  shall  henceforth  only  speak  of  stabi lizabil ity. 

It  is  important  to  find  "nice"  necessary  and  sufficient  tests 
for  stabi lizabil ity.  This  problem  appears  to  be  quite  formidable 
unless  one  specializes  to  particular  classes  (e.g.,  periodic)  of 
time-varying  systems.  We  do  have,  however,  the  following  sufficient 
condition. 

(8.9)  THEOREM.  Let  E = (F,  G,  H,  J)  over  il”(Z)  be  a linear  time- 

varying  system.  Suppose  that  one  can  find  integers  N and  K and  a 
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real  number  e > 0 with  the  following  property:  For  any  integer 

there  exists  some  with  < T < + M such  that 

(8.10)  det[R^R|(|](t^)  > e > 0 

where  R[y|  is  the  N step  reachability  matrix 

R^  :=  [G  FoG  ...  FaF  . . . a^"^Fa^“V,] 

Then,  ^ is  stabilizable. 

PROOF.  Essentially  condition  (8.10)  corresponds  to  the  system  being 
i (Z)-reachable  in  N steps  but  not  at  all  times. 

Given  any  initial  state  C in  R'^  and  any  initial  time  t^,  we 
construct  an  open-loop  stabilizing  control  law  as  follows  : we  apply 
zero  control  (i.e.  u(t)  = 0)  for  the  time  t^  < t < t^.  Then,  we 
drive  the  system  to  zero  state  in  N steps.  This  can  be  done  because 
from  (8.10)  I is  reachable  at  time  t^.  Moreover,  the  input  sequences 
applied  are  uniformly  (in  tj^)  bounded  in  norm.  Having  brought  the 
system  to  zero  state  we  apply  no  further  inputs.  We  can  thus 
stabilize  Z by  an  open-loop  control  law.  This  implies  that  Z is 
asycontrol lable,  which  by  Theorem  (8.5)  implies  that  Z is 

stabilizable.  Q 


CHAPTER  NINE 

STABLE-PROPER  FACTORIZATIONS 

Of  late,  the  use  of  stable-proper  factorizations  for  control 

system  design  has  become  increasingly  popular,  one  advantage  being 

that  properness  of  controllers  is  automatic.  See  for  example 

VIDYASAGAR  [1978],  DESOER  et  al  . [1980],  SAEKS  and  MURRAY  [1981], 

KHARGONEKAR  and  SONTAG  [1982],  etc.  In  this  chapter  we  investigate 

in  detail  stable-proper  factorizations  for  time-varying  systems. 

Throughout  this  chapter,  we  shall  treat  only  bounded  time- 

varying  systems,  i .e.  systems  defined  over  the  difference  subring 
00 

i (Z).  This  is  a reasonable  physical  constraint  as  discussed  at  the 
end  of  the  previous  chapter.  The  results  we  obtain  in  this  chapter 
extend  to  some  classes  of  time-varying  systems  defined  over  a 
difference  subring  B _ £ (Z)  such  as  periodic  systems. 

Define  the  set  RP  to  be 

\ 

RP  :=  {<1.  in  £"(Z)[[z"^]]  : = a(zl  - + 6} 

where  a,  S,  y,  and  6 are  matrices  of  compatible  dimensions  over 
00 

£ (Z).  Here,  RP  is  a mnemonic  for  _rational  or  _realizable  and 
£Toper . It  is  easy  to  verify  that  RP  forms  a ring  with  the  usual 
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00  1 

skew-multiplication  and  addition  defined  in  £ (Z)[[z~  ]].  Also, 
define  the  subring  RP^  by 

RPg  •“  + 5 in  RP  : (zl  - 3)'^  is  stable}  . 

We  shall  call  RP^  the  ring  of  stable,  proper,  rational  functions  (the 
subscript  s denotes  stable).  We  shall  abandon  rigorous  nomenclature 
and  loosely  refer  to  elements  in  RP^  or  matrices  over  RP^  as  being 

stable-proper. 

Let  D(z)  be  an  nxn  matrix  over  RP.  We  shall  say  that  D(z)  is 

bi causal  if  and  only  if  D has  an  inverse  D"^(z)  also  over  RP.  We 

have  the  following  simple  lemma. 

(9.1)  LEMMA.  Let  D(z)  = M|(zl  - M2)~^M-^  + M^  be  an  nxn  matrix  over 
RP.  Then,  D is  bicausal  if  and  only  if  M/^  is  invertible  over  t“(Z), 
and,  in  this  event  D~^(z)  is  given  by  the  formula 

(9.2)  D-l(z)  = M’^  - M’^M^(zI  - M^  + M3M'^M^) ■^M3M^  . 

PROOF.  Suppose  M4  is  invertible  over  t“(Z).  Then,  it  can  be 

verified  by  direct  multiplication  that  D"^(z)  is  given  by  (9.2).  To 

prove  the  converse,  suppose  D(z)  is  invertible  over  RP.  Let 


D"^(z)  = Nq  + N]^z‘^  + N2Z“2  + ... 
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Notice  that  DD“^  = I.  Thus, 

(M^  + terms  in  z"^,  z“^,  ...)(Nq  + N^z"^  + ...)  = I 

Equating  coefficients  of  z°  we  obtain  M^Nq  = I.  Notice  that  Nq  is 
bounded  (i.e.  over  Z (Z)).  This  proves  our  claim.  Q 

The  previous  lemma  is  merely  a generalization  of  the  familiar 
time-invariant  theory  result  that  a proper  power  series  D in  R[[z"^]] 
is  invertible  and  its  inverse  is  proper  if  and  only  if  the  constant 
coefficient  of  D is  invertible. 

A 

Let  f be  an  input/output  map  and  let  W^(z)  be  its  associated 

A 

transfer-function  matrix.  Then,  W^(z)  is  said  to  admit  a right- 
Bezout  stable-proper  factorization  if  and  only  if  there  exist  stable- 
proper  matrices  (i.e.,  over  RP^)  N,  0,  X,  and  Y with  D bicausal  and 
such  that 

(9.4)  W^(z)  = ND"^  , XN  + YD  = I 

One  can  also  similarly  define  left-Bezout  stable-proper 
factor! zati on . Such  factorizations  have  been  found  to  be  extremely 
useful  in  tackling  many  control -theoretic  problems  (for  example, 
optimal  controller  design,  see  ZAMES  and  FRANCIS  [1983]).  As  with 
the  polynomial  factorization  theory  of  Chapter  Five,  we  begin  by 
characterizing  input/output  maps  whose  transfer-functi ons  admit  left- 
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and/or  right-Bezout  stable-proper  factorizations.  We  have  the 
following  central  result. 

(9.4)  THEOREM.  Let  f be  an  input/output  map  and  let  W^(z)  be  its 
pxm  associated  transfer-function  matrix  over  &°°(Z)[[z~^]].  Then, 
the  following  are  equivalent: 

( a ) W^(z)  admits  a right-Bezout  stable-proper  factorization. 

A 

( b ) W^(z)  admits  a left-Bezout  stabl e-proper  factorization. 

(c ) f admits  a stabilizable  and  detectable  realization. 

PROOF,  (c)  > (a).  Let  Z = (F,  G,  H,  J)  be  a stabilizable  and 
detectable  realization  of  f.  Consequently,  by  definition,  there 
exist  matrices  L and  K over  £"(Z)  that  (zl  - F + GL)"^  and  (zl  - 
F + KH)“^  are  stable  matrix  power  series.  Define  the  stable-proper 
(i.e.,  over  RP^)  matrices  N,  D,  X,  and  Y by 

N = (H  - JL)(zI  - F + GL)-1g  + J , 

D = I - L(zl  - F + GL)“1g 

(9.5) 

X = L(zl  - F + KH)"k  , 


Y = I + L(zl  - F + KH)-1(G  - KJ) 
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Notice  by  that  Lemma  (9.1)  D is  bicausal.  It  can  be 
mechanically  verified  that  with  these  definitions 

W^(z)  :=  H(zl  - F)"^G  + J = ND"^  . XN  + YD  = I 

A 

Thus,  W^(z)  admits  a right-Bezout  stable-proper  factorization. 

A 

(a)  > (c).  Now  suppose  that  W^(z)  admits  a right-Bezout  stable- 
proper  factorization,  i.e.  these  exist  stable-proper  matrices  N,  D, 
X,  and  Y with  D bicausal  and  such  that 

W^(z)  = ND"^  XN  + YD  = I . 

From  the  definition  of  a stable-proper  function  (i.e.,  the  ring  RPg), 
it  follows  that  we  can  write 

N = Ni(zl  - N2)-1  N3  + N4 

D = D|^(zl  - D2)  ^Dg  + D^ 

where  ...  , N^,  D^,  ...,  D^  are  matrices  of  appropriate  sizes 
over  t (Z),  with  (zl  - N2)~^  and  (zl  - M2)~^  stable  matrix  power 
series.  Since  D is  bicausal,  it  follows  from  Lemma  (9.1)  that 
without  loss  of  generality,  D^  = I.  Define  the  linear  time-varying 
system  (over  I = (F,  G,  H,  J)  by 
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-X3°l 

‘Na' 

F = 

G = 

0 

Dj  - OjDj_ 

1 

CO 

Q 

__  1 

(9.6) 

H = [N^  - J = . 

The  system  l is  stabiiizable  because  with  L = [0  - (which  is 

over  Z°°(Z)), 


(zl  - F + GL)"^  = 


(zl  - 
0 


0 

(zl  - 


which  is  stable.  We  now  show  that  E,  defined  by  (9.6)  is  detectable. 
Define  the  stable-proper  matrix  Q by 


Q = 


(zl  - N2)'^  N3 
(zl  - D^)-^  D3 


It  is  easy  to  verify  that  (zl  - F)Q  - GD  = 0 and  HQ  + JD  = N. 
Combining  these  equations  with  XN  + YD  = I,  we  can  write 
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zI-F  G 

- _ 
(zI-F+GL)"^  Q 

I 0 

_ XH  -(Y+XJ)_ 

_L(zI-F+GL)“^  -D_ 

_V  I 

where  V is  some  stable-proper  matrix  whose  exact  formula  is  not 
critical  to  our  needs.  Define  the  matrix 


zI-F 


<l>  = 


XH 


G 

-(Y+XJ) 


It  is  clear  from  (9.7)  that  <j)  is  right-invertible  with  (the 

subscript  R denotes  right)  a stabl e-proper  matrix.  We  now  show  that 

<|)  is  also  left -invertible  (and  thus  its  inverse  is 

K L 

unique).  Notice  that 


(9.8) 


<t>  = 


zI-F 

0 


(zI-F)-^G 

-(Y+XJ) 


The  matrix  is  clearly  both  left-  and  right-  invertible.  Also, 
(Y  + XJ)  must  be  of  the  form 

(Y  + XJ)  = I + terms  in  , z"^,  ... 
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since  XN  + YD  = I.  Thus,  the  leading  (z°)  coefficient  of  (j)2  is  of 
the  form 


I 0 

* -I 


Consequently,  <p2  is  both  left-  and  right-  invertible.  Thus,  from 
(9.8)  (})  must  also  be  both  left-  and  right-  invertible,  proving  our 
claim.  Recalling  that  is  a stable-proper  matrix,  we  can  write 


* 

★ 

(zI-F)  G 

XH  -(Y+XJ) 


One  component  is  the  above  equality  is 


Wi(zl  - F)  + (W2X)FH  = I 


Thus  the  dual  of  the  pair  (F,  H)  is  asycontrollable,  which  by 
Theorem  (8.5)  implies  that  £ is  detectable.  This  completes  the  proof 
of  (a)  -»■  (c). 

(b)  (c)  and  (c)  -*•  (b)  follow  from  a dual  argument.  Q 
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(9.9)  REMARK.  In  proving  the  above  result,  we  have  made  critical  use 
of  Theorem  (8.5);  the  fundamental  result  of  the  previous  chapter. 
Recall  that  Theorem  (8.5)  states  that  asycontrol lability  is 
equivalent  to  stabilizability,  and  the  rather  difficult  proof  of  this 
result  is  in  Appendix  C.  There  appears  to  be  no  way  to  circumvent 
Theorem  (8.5)  in  studying  the  existence  of  stabl e-proper 
factorizations.  We  would  further  like  to  remark  that  equations  (9.5) 
enable  one  to  compute  stable-proper  factorizations  once  the 
stabilizing  feedback  matrices  L and  K are  determined.  Q] 

Theorem  (9.4)  essentially  states  that  a time-varying  system 
admits  stable-proper  factorizations  if  and  only  if  it  admits  a 
stabilizable  and  detectable  realization.  This  characterization  is 
non-trivial  because  (in  glaring  contrast  with  time-invariant 
systems),  not  all  time-varying  systems  admit  stabilizable 
realizations  (and  therefore  by  Theorem  (9.4),  stable-proper 
factorizations).  We  illustrate  this  with  the  following  example. 

(9.10)  EXAMPLE.  (Also  see  Example  (5.9)).  Consider  the  linear 
time-varying  system  E = (F,  G,  H)  over  4~(Z)  where  F(k)  = G(k)  = 1 
for  all  k in  Z and  G = A = the  unit  pulse  concentrated  at  the 
origin.  Let  fj.  be  the  input/output  map  associated  with  E.  Suppose 
fj;  admits  a stabl  izable  realization  E = (F,  G,  H).  This  would  imply 
that  for  any  initial  state  C at  time  t=l,  there  exists  an  open  loop 
control  law  that  drives  E from  x(l)  = C to  zero  final  state  (and 
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therefore  zero  final  output)  asymptotically.  However,  from  the  unit- 
pulse  response  function  of  fj.,  it  is  evident  that  application  of 
inputs  after  t=l  has  no  effect  on  the  output.  Thus,  it  must  be  that 

(9.11)  lim  n H(k)P(k-l)F(k-2)  ...  r(l)  n = 0 . 

k-*-« 

However, 

W.  (z)  :=  H(zl  - F)"^G  = z“^A  + z’^  A + ... 

' Z 

= fl(zl  - F)‘^G  . 


Consequently,  for  all  integers  k > 1 


R(k)F(k-l)  ...  F(1)G(0)  = 1 . 


This,  together  with  the  fact  that  5 is  over  C{1),  renders  (9.11) 
impossible.  Therefore,  fj.  does  not  admit  a stabilizable  realization 
which  implies  (by  Theorem  (9.4))  that  W,  (z)  does  not  admit  a 


stable-proper  factorization. 


CHAPTER  TEN 
FEEDBACK  CONTROL 

Having  investigated  in  detail  issues  related  to  the  existence  of 
stable-proper  factorizations  in  the  previous  chapter,  we  now  examine 
their  role  in  the  problem  of  stabilization. 

Let  I = (F,  G,  H,  J)  and  = (F^.,  G^,  H^,  J^)  be  linear  time- 
varying  systems  over  a"(Z)  called  the  plant  and  controller 
respectively.  Consider  the  feedback  configuration  shown  below. 


v(k) 


■>y(k) 


We  have  the  following  result. 

(10.1)  THEOREM.  The  controller  internally  stabilizes  z if  and 
only  if  the  following  conditions  are  met: 

(a)  Both  Z and  Z^  are  stabilizable  and  detectable. 

(b)  Let 
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Wj.(z)  :=  H{zl  - F)"^G  = D'^N  , NX  + DY  = I 


(z)  :=  H (zl  - F )"^G  =ND'^,XN+YD=I 
c'  c'  c cc  cc  cc 

c 


be  any  stable-proper  factorizations  of  the  plant  and 
controller  transfer-function  matrices  (these  exist  by 
(a)).  Then, 


S”1  ;=  (DOj.  + NN^)-1 
is  a stable  matrix  power  series. 

SKETCH  OF  PROOF:  The  proof  of  the  above  theorem  is  tedious,  but 

otherwise  straightforward.  Observe  first  that  internally 

stabilizes  E if  and  only  if 

(10.2)  : = 

is  a stable  power  series. 

Suppose  E^  stabilizes  E.  Then,  (a)  follows  directly  fron 
Corollary  (8.7).  To  prove  (b)  one  first  computes  Each  of  the 

four  block  submatrices  of  (}>"^  is  stable  and  these  four  conditions 
imply  (after  a series  of  manipulations)  that  S“^  is  stable. 


zI-F  -GH 


zI-F 
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Now  suppose  conditions  (a)  and  (b)  are  met.  To  prove  that 
internally  stablizes  I one  mechanically  verifies  that  with 
stable,  each  of  the  four  (block)  submatrices  of  <|)~^  are  also  stable. 

(10.3)  REMARK.  Essentially,  the  above  theorem  permits  us  to  design 
controllers  (that  i nternal ly  stabilize  the  plant)  working  exclusively 
in  the  "frequency-domain".  Internal  stability  is  realization 
dependent,  and.  Theorem  (10.1)  says  in  effect  that  as  long  as  we 
i mpl ement  stabilizable  and  detectable  realizations  of  our 
controllers,  internal  stability  of  the  closed-loop  system  is 
guaranteed  if  (and  only  if)  S“^  is  a stable  matrix  power  series.  D 

We  can  now  parameterize  all  controllers  that  stablize  a given 
plant  2.  We  shall  assume  that  D is  stabilizable  and  detectable  (and 
ther'efore  admits  stable-proper  factorizations);  else  by  Theorem 
(10.1)  Z cannot  be  stabilized.  We  have  the  following  theorem  (the 
proof  closely  resembles  that  of  DESOER  et  al . [1980,  Theorem  3]  and 
is  therefore  omitted). 

(10.4)  THEOREM.  Let  Z = (F,  G,  H,  J)  be  a stabilizable  and 
detectable  plant.  Let 

W^(z)  = ND'^  = D“^Nj  , 

XN  + YD  = I , + D^Yj  = I 
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be  any  left-  and  right-  stable-proper  factorizations  of  W^.  Then,  a 

controller  internally  stabilizes  £ if  and  only  if  the  controller 

transfer-function  Wy  (z)  is  of  the  form 

^c 

(10.5)  (z)  = (MN^  + Y)"^  (-MD^  + X) 

for  some  matrix  M over  RP^ . Q] 

Note  the  close  resemblance  (10.5)  bears  to  the  time-invariant 
case  (DESOER  et  al  . [1980,  Theorem  3]). 

Armed  with  the  powerful  tool  of  stable-proper  factorizations,  we 
take  a brief  look  at  track!  ng-probl  ems  for  linear  time-varying 
systems . 

Consider  the  feedback  system  configuration  shown  below. 


Here,  Z is  the  (given)  plant,  is  the  (to  be  designed)  controller, 
and  Z[^  is  the  (given)  exogeneous  system  that  generates  the  reference 
signal  r(k)  via  initial  conditions  Xp(kg)  on  Z|^.  The  system  Zp  could 
in  general  be  time-varying.  For  example,  Z[^  could  generate  a 
reference  signal  r(k)  which  is  a sinusoid  of  time-dependent 
frequency,  or  r(k)  could  be  a polynomial  in  time. 
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(10.6)  EXAMPLE.  Suppose  = (Fp,  Gp,  Hp)  is  described  by 


F,(k) 


cos(TQ+Tj^+2kT^) 

sin(TQ+T^+2kT^) 


-sin(TQ+T^+2kT^)' 
cos (TQ+T^+2kT^) 


HpCk)  = [1  0] 


Gr(K) 


where  Tg  and  are  (fixed)  real  numbers.  Then,  for  the  initial 
condition  x(0)  = [A  0]',  the  reference  signal  r(k)  generated  by 
Sr  is 


r(k) 


0 k<0 

Acos(kT(k)),  k>0 


where  T(k)  = Tg  + kT^.  Thus,  in  this  case  the  reference  signal  r(k) 
is  a sinusoid  whose  frequency  is  modulated  by  a ramp.  Q 


We  consider  the  following  design  problem,  called  Tracking  with 


Internal  Stability  (TIS): 
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Find  (if  possible)  a controller  such  that 

(i)  internally  stabilizes  E. 

(ii)  For  any  initial  time  kg  and  any  initial  state  C,  the 
output  y(k)  of  the  closed  loop  system  asymptotically 
tracks  the  reference  signal  r(k)  generated  by  Ej^  with 
initial  condition  Xp(kg)  = 5)  uniformly  in  Kg. 

The  TIS  problem  can  be  reformulated  in  the  "frequency  domain"  as 
follows.  Let  P(z),  C(z),  and  T(z)  denote  the  plant,  controller,  and 
exogenous  system  transfer  function.  Then,  find  (if  possible)  C(z) 


that 

(i) 

(I 

+ PC)“^  is  a stable  matrix  power  series. 

(ii) 

(I 

+ PC)"^T  is  a stable  matrix  power  series. 

Observe  that  by  Theorem  (10.1)  E^  stabilizes  E internally  only 
if  both  E and  E^  are  stabilizable  and  detectable.  Hence,  by 
Theorem  (9.4),  the  plant  and  controller  tranfer-functions  P(z)  and 
C(z)  must  admit  stable-proper  factorizations  in  order  to  meet 
requirement  (i).  We  further  assume  that  the  exogenous  system 
transfer-function  T(z)  admits  a left-Bezout  stabl e-proper 
factorization.  In  the  event  that  Ej^  is  time-invariant,  obviously 
such  a factorization  always  exists. 
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SAEKS  and  MURRAY  [1981]  have  formulated  the  TIS  problem  in  a 
general  axiomatic  setting  based  on  stable-proper  factorizations.  The 
TIS  problem  for  time -varying  systems  falls  into  this  general 
framework.  The  difficulty  involved  however,  is  in  showing  the 
existence  of  stable-proper  factorizations  for  time-varying  systems 
and  this  was  done  in  the  previous  chapter.  In  terms  of  these 
factorizations  following  SAEKS  and  MURRAY  [1981,  Theorem  2]  we  give  a 
partial  solution  to  the  TIS  problem. 

Let  us  consider  left-  and  right-  stable-proper  facotrizations  of 
the  plant  and  exogenous  system  transfer-functions  as 

P(z)  = NO"^  = 

(10.7)  XN  + YD  = I = + D^Y^ 

T(z)  = W-^V  , Wy^  + VU2  = I . 

Then,  we  have  the  following  result  (which  follows  directly  from  SAEKS 
and  MURRAY  [1981,  Theorem  2]. 

(10.8)  THEOREM.  The  problem  TIS  admits  a solution  if  and  only  if 
there  exist  stable-proper  matrices  M^  and  M2$uch  that 

(10.9)  NMjDj  + M2W  = NX  - I 

In  this  event,  the  set  of  all  controllers  that  solve  the  TIS  problem 
can  be  parameterized  as 
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C(z)  = (M^Nj  + Y)-l(-MiDi  + X) 

where  is  any  stable-proper  solution  of  (10.9).  | [ 

It  is  desirable  to  have  nice  necessary  and  sufficient  conditions 
for  the  solvability  of  the  linear  skew  equation  (10.9),  and,  in  the 
event  that  (10.9)  is  solvable,  one  would  desire  constructive 
procedures  to  compute  the  solution.  This  open  problem  appears  to  be 
quite  formidable.  We  would  also  like  to  remark  in  closing  that,  as 
is  well  known,  the  problem  of  disturbance  rejection  can  be 
reformulated  as  a TIS  problem  and  handled  similarly. 


CHAPTER  ELEVEN 
CONCLUDING  REMARKS 

In  this  dissertation  we  have  developed  a "transfer-function" 
type  theory  for  linear  time-varying  discrete-time  systems.  Using 
this  framework,  in  the  first  part  of  the  dissertation,  we  have  been 
able  to  generalize  much  of  the  existing  polynomial  model  theory. 
Specifically  we  have  treated  polynomial  realization  theory 
(FUHRMANN),  polynomial  factorization  theory,  and  applications  to 
feedback  control.  In  the  second  half  of  the  dissertation,  we  have 
treated  the  problems  of  stabilization,  of  existence  of  stable-proper 
factorizations,  and  taken  a cursory  look  at  the  tracking  problem  for 
time-varying  systems.  One  of  our  deepest  results  is  the  equivalence 
of  dynamic  and  memoryless  state  feedback  as  far  as  the  problem  of 
stabilization  is  concerned. 

Many  open  problems  remain  in  the  area  of  time-varying  systems. 
We  outline  some  of  these  below,  and  expect  that  the  tranfer-function 
theory  developed  in  this  dissertation  will  prove  useful  in  solving 
some  of  these  problems. 

(a)  Regulation  and  Tracking.  We  have  shown  in  Chapter  10  that 
the  problem  of  Tracking  with  Internal  Stability  reduces  to 
the  solution  of  a linear  matrix  equation  (10.9)  over  the 
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skew-ring  RP^.  It  would  be  very  desirable  to  give 
explicit  techniques  for  the  solution  of  such  equations, 
perhaps  for  a more  restricted  class  of  time-varying 
systems.  It  would  be  also  interesting  if  one  can  extract 
a form  of  the  Internal  Model  Principle  (WONHAM  [1979])  in 
this  framework. 

(b)  Many  of  the  results  of  this  dissertation  are  not  concrete 
enough  that  they  may  be  called  "completely  computable". 
This  is  because  we  treat  completely  general  time-varying 
systems.  To  use  the  theory  to  design  and  implement 
controllers  one  has  to  obtain  much  more  computable 
results,  and,  perhaps  an  extended  study  of  particular 
classes  of  time-varying  systems  (such  as  periodic,  or 
slowly-varying)  may  yield  such  results.  Specifically,  we 
would  like  to  solve  the  sensitivity  minization  problem  of 
ZAMES  and  FRANCIS  [1983]  and  any  of  the  numerous  problems 
in  robust  multivariable  control  system  design  (see  DOYLE 
and  STEIN  [1981]  for  example)  for  these  particular  classes 
of  time-varying  systems. 

(c)  ANDERSON  and  MOORE  [1981b]  have  pointed  out  that  for  time- 
invariant  large  scale  systems  it  is  possible  to  eliminate 
"fixed -modes"  by  using  decentralized  time-varyi ng 
feedback.  It  thus  appears  that  the  use  of  time-varying 
controllers  for  time-invariant  plants  may  be  advantageous. 
A detailed  study  of  this  problem  would  prove  quite 
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interesting.  For  example  can  one  improve  stability 

margins  by  the  use  of  time-varying  feedback,  what  remains 
invariant  under  alj_  types  of  feedback,  etc.? 

(d)  As  an  extended  project,  one  might  consider  continuous-time 
generalizations  to  the  transfer-function  theory  developed 
in  this  dissertation.  This  problem  is  quite  non-trivial 
as  much  of  the  dissertation  inextricably  uses  the 
"discrete-time  nature"  of  the  problem.  For  instance  the 
action  of  z is  a right  time-shift  of  the  coefficients  for 
discrete-time  systems,  and  is  a derivation  for  continuous 
time  systems.  The  interested  reader  may  consult  KAMEN 

[1976]  for  related  results  on  continuous-time  time-varying 
systems. 


APPENDIX  A 

PROOF  OF  PROPOSITION  (4.1) 

Before  we  establish  Proposition  (4.1),  we  shall  require  several 

n 

intermediate  lemmas.  Let  Q(z)  = be  an  r x r invertible 

polynomial  matrix  over  A[z]  and  let  its  inverse  over  A((z"^))  be 

Q ^(z)  = j=^N  ^ ^ Qi  •=  0.  Define  a right 

A-1  inear  projection  map  ttq  (as  before)  by 

Tfq  : a’'[z]  > A^[z]  : x > Qtt{Q'^x} 

Also,  define  a left  A-1 inear  projection  map  Xg  by 

: A^[z]  A^[z]  : y -»■  ir{yQ‘bQ 

Here,  elements  y in  A^[z]  are  to  be  thought  of  as  row  vectors. 
Notice  that  TTg(x)  and  Xg(y)  are  polynomial  vectors  of  degree  less 
than  n.  We  can  therefore  write 

(A.l)  TT  (iz*^)  = E z^  <|)  , 

H j=0 
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4 


Xp(lz') 


n-1 

Z 

j=0 


Notice  that 


7Tq(Iz^)  = Qtt{Q-^z''} 


n 

= E 
i=0 


1=1 


-i 


n-1 

E 

j=0 


Equating  powers  of  z,  we  get  after  some  computation. 


(A. 2) 


n 

= E 

t=l 


8t.k> 


Similarly,  we  can  obtain 


(A. 3) 


♦j.k 


A A 


Now  define  the  nr  x nr  block  matrices  Q,  B,  M,  and  N over  A by 


(A. 4) 


l,j-l 


N.  . :=  ip.  , . , 
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(A. 5)  LEMMA.  With  the  above  definitions  for  Q,  B,  M,  and  N, 
M = Q B . N = a(B  Q) 

PROOF.  We  verify  the  first  equation  by  direct  computation: 


■ tfi  '^tj  • 


Thus,  M = Q B.  The  second  equation  can  also  be  easily  verified.  Q 


Notice  now  that 


Q(Q'^z^)^  = I - 7T  (Iz*") 


is  a polynomial  matrix  of  degree  < n.  Thus,  we  can  write 


j=0 


n-1 


(A. 6) 


(z'^Q'^)  Q = E 0.  .zJ. 
j=0 


n-1 
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A direct  computation  for  a and  g (in  a manner  similar  to  obtaining 
(A. 2),  (A. 3))  yields 


X • 1 

J,k 


n 

= E 
t=l 


a" 


-t 


^^t+j-n  \+k-n- 


(A. 7) 


n 

E 

t=l 


(B 


t+k-n 


) 


(Q 


t+j-n 


) 


Define  the  nr  x nr  block  matrices  X,  Y,  U,  and  V over  A by 


^^i+j-n-1^’ 


Y.  . 

i.J 


(B...  J 
^ i+j-n-r 


(A. 8) 


a. 


V.  . 
1 ,J 


= B 


(A. 9)  LEMMA.  With  the  above  definitions. 


U = XY  , a"’\v)  = YX  , 

M+U=I  , N+V=I,  M^=M 


no 


PROOF.  The  first  two  equations  can  easily  be  verified  by 
computation.  We  derive  the  third  equation: 

Notice  that 


7r^(z''l)  + (Q  Q-1  z'')^  = z^I 


n-1 
= z 
j=0 


* “j,k> 


Thus, 


j=k 

otherwi se 


It  now  follows  from  the  definitions  of  M and  U and  M + U = 
identical  fashion,  we  can  obtain  N + V = I. 

We  can  now  prove  that  = M.  Notice  that 


n-1 


n-1 


a di rect 


I.  In  an 


The  right  A-1  inear  map  ttq  is  a projection.  Consequently, 


Ill 


j=0  t=0  ^t+l,j+l  ^>l,k+l 


Equating  powers  of  z,  we  get 

n-1 

^+l,k+l  " ^+l.s+l  ^s+l,k+l* 

i.e.,  = M.  This  completes  the  proof.  D 

We  shall  also  require  the 

(A. 10)  PROPOSITION.  Xq  is  a finitely  generated  right  A-modul e. 

Further,  let  M be  the  nr  x nr  matrix  described  above.  Then  Xg  is 

isomorphic  to  Image  (M)  (as  a right  A-modul e) . 

PROOF.  We  first  prove  that  Xg  is  finitely  generated.  Since  Q is  an 
r X r invertible  polynomial  matrix,  it  follows  from  Proposition  (3.2) 
that  there  exists  an  r x r polynomial  matrix  T such  that  QT=M  is 
monic  and  such  that  deg  (Q)  = deg(M)  = n.  For  any  x in  Xg, 
Q'^x  =ir(Q"lx).  Therefore  deg  (T^^Q-^x)  < deg  (T“^)-l.  Notice  that 
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T ^ is  M ^Q.  Thus,  deg  (T"^)  = deg(Q)  - deg(M)  since  M is  monic. 
However,  deg(Q)  = deg(M),  therefore,  deg  (T-^)  = 0,  i.e.,  is 
proper.  Consequently, 

-n  + deg(x)  = deg  (M“^x)  < deg  (T"^)  - l = _i 

Thus,  for  any  x in  Xq,  deg(x)  < n-1.  Thus,  Xg  is  finitely  generated 
and 


‘N.J  ^ '"’-I- 


n-1 


1,  2, 


r} 


where  ej  is  the  j-th  column  of  the  r x r identity  matrix,  generate 
Xq. 

We  now  show  that  Xg  = Image  (M).  Consider  the  map 


n-1  r 


<|)  : X^  + Image  (M)  : E E b.  . a.  .->■  M a 

i=0  j=l 

where  ^ = [ag^^  “0,rl  “1,1  “l,rl  •••  ]“n-l,l  ••• 

^n-l,r^  ’ straightforward  to  verify  that  (|)  is  well-defined  and 

right  A-1  inear. 

Also,  since  generate  Xg,  <j)  is  onto.  Suppose  that 

<P(x)  = 0,  for  some  x in  Xg.  Let  ^ be  partitioned  as  a 

= [aj  ...  a'  i]‘  where  each  a.  is  an  r x 1 block.  Then,  Ma  = 0, 
n-1  . — 


i.e.,  TTp.  ^ E z a*)  — 0. 
^ i =0  ^ 
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Notice  however  that 


n-1  r 


n-1 


X = E E b.  . a.  . = E 

i=0  j=l  ■'’J  i=o 

Thus,  Ip  is  one-to-one.  Consequently,  Xg 
This  completes  the  proof. 


Tfg(z'')  a.  = 0 
= Image  (M). 


□ 


The  following  lemma  is  a simple  (discrete-time)  version  of  a 
theorem  of  Dolezal  [1964]. 

(A.  11)  LEMMA.  Let  M be  any  p x q matrix  over  A.  Then,  the  right 
A-module  Image  (M)  is  free  if  and  only  if  rank  M(t)  = p = constant 
for  all  t in  Z. 

PROOF.  The  proof,  being  straightforward,  is  omitted.  D 


We  are  now  in  a position  to  prove  the 
(4.1)  PROPOSITION  Xn  is  a finitely  generated,  free  right  A-module. 
PROOF.  We  have  already  shown  that  Xg  is  finitely  generated  (Lemma 
(A.IO)).  It  follows  from  Lemmas  (A. 10),  (A.l)  that  in  order  to  prove 
Xg  is  free,  we  need  only  to  show  that  rank  M(t)  = p = constant,  for 
all  t in  Z.  We  proceed  to  do  this. 

Notice  that  from  Lemma  (A. 5) 


det  [XI  - M(t)]  = det  [XI  - Q(t)B(t)] 
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From  Lemma  (A. 9),  we  see  that  M = I - U = I - XY.  Consequently, 
det  [XI  - M(t)]  = det  [(X-l)I  + X(t)Y(t)] 

= det  [(X-l)I  + Y(t)  X(t)  . 

Again  from  Lemma  (A. 9),  we  have  YX  - I = a^'^(V-I)  = - a^“l(N). 
Therefore, 

det  [XI  - M(t)]  = det  [XI  - N(t  - n+1)] 

It  now  follows  from  Lemma  (A. 5)  that 

det[XI  - M(t)]  = det  [XI  - B(t-n)Q(t-n)] 

= det  [XI  - Q(t-n)B(t-n)] 

= det  [XI  - M(t-n)] 

Therefore,  the  eigenvalues  of  M(t)  are  identical  to  the  eigenvalues 
of  M(t-n).  Since  M(t)  is  idempotent,  i.e.,  w2(t)  = M(t)  (see  lemma 
(A. 9)),  it  follows  from  standard  linear  algebra  that  M(t)  is  simple. 
Thus,  rank  M(t)  = number  of  eigenvalues  that  are  zero,  which  in  turn 
implies  that 


(A. 12) 


rank  M(t)  = rank  M(t-n) 
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We  could  also  write 

Q(z)  = + H z‘'o.  = 0*(z) 

i=0  ^ 

Corresponding  to  Q*,  we  obtain  an  (n+1)  r x (n+l)r  matrix  M*  in  a 
manner  analogous  to  the  way  M was  obtained  from  Q.  Notice  that  by 
Lemma  (A. 10) 


Xq  = Image  (M)  = Image  (M*). 

Thus,  for  all  t in  Z,  rank  M(t)  = rank  M*(t).  By  an  argument 

identical  to  the  one  presented  above,  we  can  conclude  that  for  all  t 
in  Z 


rank  M*(t)  = rank  M*(t-n-l)  , 

i.e.,  rank  M(t)  = rank  M(t-n-l).  This  together  with  (A. 12)  implies 
that  rank  M(t)  = p = constant  for  all  t in  Z.  I I 


APPENDIX  B 

PROOFS  OF  THEOREMS  (5.2)  AND  (5.11) 

We  would  like  to  first  emphasize  that  the  proofs  of  Theorems 
(5.2)  and  (5.11)  given  in  the  following  pages  is  constructi ve.  A 
systematic  procedure  for  obtaining  the  polynomial  factorizations  of 
these  theorems  is  given  in  Chapter  5,  and  essentially  follows  from 
this  appendix.  Before  we  prove  these  theorems,  we  need  to  establish 
a few  intermediate  results,  some  of  which  are  themselves  of  interest. 
We  proceed  to  do  this. 

Let  (F,G)  be  a pair  of  matrices  over  A that  are  reachable  in  N 
steps  at  all  times.  It  follows  from  Lemma  (2.4)  that  the  matrix 

R^  :=  [G|F(aG)|...|F(oF)(a2F)  ...  (o^"^F) (a^'^G) ] 

is  right-invertible  over  A.  This  implies  that  there  exist  matrices  A 
and  C such  that 


[-F  G] 


C 
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i.e.,  the  row  [-F  G]  is  unimodular  (over  A).  It  is  now  a classical 
fact  that  there  exist  matrices  X and  Y such  that 

'-F  g" 

M = 

0 

X Y 

is  unimodular  (see  NEWMAN  [1973.  application  (d),  p.  38]). 

We  now  have  the  following  key  result. 

(B.l)  PROPOSITION.  Let  (F,  G)  be  a pair  of  matrices  over  A that  are 
reachable  in  N steps  at  all  times,  for  some  integer  N.  Let 


-F  G 

M = 

0 

_ X Y_ 

be  any  unimodular  completion  of  the  row  [-F  G]  over  A and  let  N^  : = 
Mq~^  be  partitioned  as 


A b' 

N = 

0 

C 0 

where  A is  an  nxn  matrix.  Then,  the  pair  (A' , B* ) is  observable  in  K 
steps  for  some  integer  K. 
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PROOF.  Since  the  pair  (F,  G)  is  reachable  in  N steps  at  all  times, 
it  follows  from  Proposition  (2.8)  that  there  exists  an  m x n matrix  L 

A 

over  A such  that  F :=  F + GL  is  a-nil potent,  i.e.,  for  some 
integer  k 

F(aF)  ...  (a^  F)  = 0. 


Notice  that 


Here  the  precise  formulae  for  C,  D,  and  X are  not  critical  to  our 
needs.  One  component  in  the  above  equality  is 


(B.2) 


- AF  + BX  = I 


We  now  see  that 
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F = (a-^  I)  F 

= - (a’^A)(a’^  F)  F + (o'^B)(a“^  X)  F . 
Substituting  this  into  (B.2)  yields 

(B.3)  BX  - A(a'^  B)(a"^  X)  F + A(a‘^A) (a“^F)  F = I 

From  (B.2)  it  follows  that 
(a'^  F)  F = (a‘^  I)(a'^  F)  F 

= - (a'2  A)(a'^  F)(a‘^  F)F  + B)(a'^  X)(a‘^  F) 

Using  this  in  (B.3)  gives  us 

BX  - A(o"^  B)(a’^  X)F  + A(a’^A)(a~^  B)(a"^  X)(a’^  F)  F 

- ^{a~^  A)(a"^  A)(o‘^  F)(a‘^  F)  F = I. 


Repeating  this  procedure  k times,  we  obtain 
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[B|A(a"^B)|  ...|A(a'^A)  ..  ^A)(a'^  B)] 


\-U 


-k  ; 


+ S,  (a  F)  ...  (a'"F)F  = I , 


where  Sg,  S^,  ...,  S|^  are  matrices  over  A whose  exact  formulae  are 
not  critical  to  our  needs.  Notice  now  that 


...  (a"-^  F)  F = 0 


because  F is  0-nil  potent.  Hence  the  matrix 


[B|  A(0"^B)|  ...|A(a"^A)  ..  (0"‘^'^^A)  (a'^B)  ] 


is  right-invertible  over  A.  It  is  then  apparent  (see  WEISS  [1972] 
and  also  Chapter  2)  that  the  pair  (A',  B')  is  observable  in  k steps 
at  all  times. 


As  an  immediate  consequence  of  this  proposition,  we  have  the 
(B.4)  LEMMA.  (F,  G)  be  a pair  of  matrices  that  are  reachable  in 

N steps . Let  A and  B be  matrices  as  defined  in  Proposition  (B.l). 
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Then,  there  exists  an  mxn  matrix  S such  that  A r=  A + BS  is 
g~^-ni1  potent,  i.e.,  for  some  integer  k 

A (a~^  A)  ...  A)  = 0. 


PROOF.  Define  the  matrices  A^  and  over  A by 
\{t)  = A(-t)  , B^(t)  = B(-t). 

It  follows  from  Proposition  (B.l)  that  the  pair  (A‘,  B')  is 

observable  in  N steps,  for  some  integer  N.  Hence,  by  duality,  the 
pair  (Ap  B^^)  is  reachable  in  N steps.  We  now  apply  Proposition 
(2.8)  and  conclude  that  there  exists  an  mxn  matrix  S-^  over  A such 

A 

that  A^  :=  A^  is  o-nilpotent,  i.e.,  there  exists  an  integer  k 

such  that  for  any  time  t 

Aj(-t)  A^(-t  - 1)  ...  A^(-t  - k)  = 0. 

Define  a new  matrix  S over  A by  S(t)  = S^C-t).  Then  for  any  time  t. 


A(t)  A(t  + 1)  ...  A(t  + k)  = 0, 


□ 


where  A :=  A + BS. 
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(B.5)  LEMMA.  Let  (F,  G)  be  a pair  of  matrices  that  are  reachable  in 
N steps  at  all  times.  Then,  there  exist  polynomial  matrices  P(z)  and 


Q(z)  A[z]  with  0(z)  monic  and  such  that 


(zl  - F)“^  G = P(z)  Q-\z). 

PROOF.  Since  the  pair  (F,  G)  is  reachable  in  N steps  at  all  times, 
it  follows  that  there  exist  matrices  Aq,  A^,  ...,  over  A such 

that 

N-1 

2 M.  A.  + I = 0. 
i =0  ^ ^ 

where  the  matrices  M^  are  defined  recursivly  by  Mq  = G,  M^-^^  = 

F(oM^).  Multiplying  the  above  equation  on  the  right  by  M^j  yields 
/ 

N 

E M.  B.  = 0 
i =0  ^ 


where  B^-  :=  A^-  M^j  for  i =0,  1,  ...,  N - 1,  and  B,sj  :=  I.  Define  a 
monic  polynomial  matrix 


N 

Q(z)  :=  E B. 

i=0  ^ 


It  is  an  easy  computation  to  verify  that  Tr{(zl  - F)"^  G Q}  = 0. 

Consequently,  (zl  -F)"^  G Q = P is  a polynomial  matrix,  verifying  our 
claim.  Q] 
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We  are  now  in  a position  to  prove  the  following  result. 

(5.11)  THEOREM.  Let  Z = (F,  G,  H)  be  a linear  time-varying  system 
that  is  reachable  in  N steps  at  all  times.  Then,  there  exist 
polynomial  matrices  P and  Q over  A[z]  where  Q is  invertible  such  that 

Wj.(z)  = H(zl  - F)"^  G = PQ"^ 

and  the  Fuhrmann  realization  E(P,  Q,  I)  is  A-isomorphic  to  E = 
(F,  G,  H). 

PROOF.  In  order  to  prove  this  result,  we  have  to  find  matrices  C,  D, 
Yj  ...  Yg  such  that  the  isomorphism  conditions  of  Theorem  (4.20)  are 
satisfied. 

Since  the  pair  (F,  G)  is  reachable  in  N steps  at  all  times,  the 
unimodular  row  [-F  G]  can  be  completed  to  a unimodular  (invertible 
over  A)  matrix  (see  discussion  before  Proposition  (B.l)).  Let 


be  any  such  completion,  and  let 
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Nq  :=  Mq 


B 

D 


where  A is  an  n x n matrix.  From  Lemma  (B.4)  it  follows  that  there 
exists  an  m x n matrix  S over  A such  that  A :=  A + BS  is 
o ^-nilpotent,  i.e.,  for  some  integer  k. 


A o‘^(A)  ...  a-^^l(A)  = 0. 


Define  the  matrix 


(aS) 


0 

0 . 


We  claim  that  the  polynomial  matrix  U :=  zW  + Mq  is  unimodular 

k-i 

(invertible)  over  A[z]  and  its  inverse  is  V :=  Nq  ( Z (-zWN„)''].  We 

. - . i =0  ^ 

verify  this  by  direct  computation.  First,  notice  that 


" I 0 

A B 

A B 

(oS)  0 

C D 

(aS)A  (aS)B  _ 

1 

o 

1 

r -n 

A 8 

_{aS)  0 

^0  0 
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Therefore 


(zWNg)'' 


'i  o' 

{ 

z 

> 

iaS)  0 

0 0 

1 

k-l|- 

I 0 

A 

B 

z 

I 0 

Q 

00 

o 

1 

0 

0 

{aS)  0_ 

k-1 


’i  o' 

A B 

’ I o' 

” 

z 

1 

z 

(aS)  0 

1 

0 0 

S 0 

- 

= z 


' I 

o' 

(aS) 

0 

(Az) 

0 


k-1 


Observe  that 
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(Az)^'^  = A (a'^A)  ...  (a'^'^^A) 


because  A is  a"^  -niipotent. 

= 0 , 

Consequently,  (-zWNq)*^  = 0.  It  then  follows  that 


k-1 

UV  = (zW  + M ) N.  ( E (-zWNj') 
^ ^ i =0  ^ 


k-1 

= (zWN.  + I)  ( E (-zWN„)M 
i =0 

= I + (-zUNq)*"  = I , 

proving  our  claim. 

We  now  partition  V as 


D 

-Q 


where  Q is  an  mxm  polynomial  matrix  over  A[z].  Let  us 
P :=  HD.  Let  Y3  = z(aS)  + X.  Then 


(B.6) 


D 


= I 


-Q 


def i ne 
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Since  the  pair  (F,  G)  is  reachable  in  N steps  at  all  times,  it 
follows  from  Lemma  (B.5)  that  there  exist  polynomial  matrices  P 
and  Q where  Q is  monic  such  that  (zl  - F)“^  G = P Q"^.  Therefore, 


where  Yy  is  some  polynomial  matrix.  We  now  see  that 


r ^ n 

— — — 

r “ - 

0 

p 

p 

= vu 

= I 

.'7. 

_-Q_ 

_-o  _ 

Hence, 


— 1 
o 

1 

II 

1 

-< 

ro 

r 

o 

I 

o 

1 

1 

1 Q. 

1 

>- 

1 

>- 

1 

O' 

1 

r 

-< 

L- 

n 

O 1 
1 

Consequently,  QY^  = Q,  which  from  Proposition  (3.2)  allows  us  to 
conclude  that  Q is  right-invertible.  Let  ip  in  ((z"^))  be  any 

right-inverse  of  Q,  i.e.,  Q»|^  = I.  To  see  that  Q is  also  left- 
invertible,  we  examine  (B.6).  Two  components  of  this  equation  are 
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(B.7) 

(zl  - F)  D = GQ 

(B.8) 

Y3D  + Y4Q  = I 

Multiplying  (B.7)  by  Y3(zl  - F)"^  on  the  left  we  obtain 
Y3D  = Y3(zI  - F)"^  GQ  . 

Substituting  this  into  (B.8)  we  see  that 
I = (Y3(zI  - F)-^G  + Y4)  0 , 

i.e.,  Q is  also  left-invertible.  Hence,  Q has  a unique  inverse  Q“^ 
in  A'^x^((z-1)). 

We  have  thus  constructed  polynomial  matrices  P and  Q where  Q is 
invertible,  such  that 

P = HD  , PQ-1  = H(zl  - F)"1g  , 

(B.9) 

GYj  + (zl  - F)Y2  = I , Y3D  + Y4Q  = I. 

It  is  now  straightforward  to  check  that  the  isomorphism  conditions  of 
Theorem  (4.20)  are  satisfied  by  choosing  C = G,  Yg  = Yg  = 0 and  with 
Yi»  Y2,  Y3,  Y^,  and  D as  in  (B.9).  This  completes  the  proof.  Q 
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(B.IO)  REMARK.  If  the  system  E (F,  G,  H)  is  assumed  to  be  canonical 

instead  of  only  reachable  in  N steps  at  all  times,  the  isomorphism 
conditions  of  Theorem  (4.20)  reduce  to  (much  simpler  form)  the 
existence  of  polynomial  matrices  Y3,  such  that 

(B.ll)  Y3P  + Y4Q  = I , H(zl  - F)-1q  = PQ-1  . 


Thus,  Theorem  (5.11)  tells  us  that  if  E = (F,  G,  H)  is  a canonical 
linear  time -varying  system,  there  exist  polynomial  matrices  P and  Q 
where  Q is  invertible  such  that  (B.ll)  is  satisfied  for  some  Y3,  Y^.n 
In  view  of  the  above  remark,  we  have  the 


(5.2)  THEOREM.  Let  f be  an  input-output  map  and  let  W^(z) 


= E A.  z"-^ 
j=l  ' 


be  the  pxm  transfer  matrix  associated  with  f.  Then  the  following 
statements  are  equivalent. 

(a)  There  exist  polynomial  matrices  P,  Q,  Y^,  an^  Y2  where  Q 
is  invertible  such  that 


W^(z)  = PQ'^  , Y^P  + Y^Q  = I 

A 

i . e . , ( z ) admits  a left-Bezout  polynomial  factorization, 

(b)  There  exist  polynomial  matrices  Q^,  R,  Y3,  an^  Y4  where  Q 
is  invertible  such  that 


W^(z)  = , Q^Y3  + RY^  = 


I 
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i.e.,  W^(z)  admits  a right-Bezout  polynomial  factoriza- 
tion . 

( c ) f admits  a canonical  realization. 

SKETCH  OF  PROOF:  (c)  (a)  is  a restatement  of  Remark  (B.IO).  (a)  -*■ 

(c).  Let  i:(P,  Q,  I)  be  the  Fuhrmann  realization  associated  with 

polynomial  representation  W^(z)  = PQ"^.  From  Theorem  (4.12)  since 
R - I,  it  follows  that  Z is  reachable  in  N steps,  and  from 
Remark  (4.13),  it  follows  that  E is  observable  in  N steps  (because 
YjP  + Y2Q  = I).  Therefore,  E is  a canonical  realization  of  f. 

(a)  -*•  (b)  follows  from  duality.  I 1 


APPENDIX  C 

PROOF  OF  THEOREM  (8.5) 


Before  we  begin  proving  Theorem  (8.5),  we  first  give  a precise 
system-theoretic  interpretation  of  asycontrollability  which  can  be 
regarded  as  a definition  in  lieu  of  Definition  (8.2).  The 
equivalence  of  these  definitions  is  demonstrated  in  KHARGONEKAR  and 
POOLLA  [1984B]. 


C.l  DEFINITION.  Let  e > 0 be  any  (fixed)  real  number.  A system 
^ = (F,  G,  H)  over  t°°(Z)  is  said  to  be  asycontrollable  if  and  only  if 

A 

there  exists  a real  number  M and  an  integer  N such  that  for  any 

initial  time  tn  in  Z and  any  initial  state  g in  R'^  with  ogn  = 1, 

there  exists  an  input  sequence  u.  (to),  u.  (to+1),  ..., 

^tn  which  results  in  the  state  trajectory  (to)  = C, 

U,^  — u 

^to  •••’  ^to  and  such  that 

0,5 


Dx.  (tp,+N)H  < e 
^0.5  ^ 


Hx  (t.+k)n  , Hu  (t„+k)ii  < M . □ 

0.5  ’^0.5  ° 
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The  subscripts  tQ  and  5 on  the  inputs  u and  states  x signify 
their  dependence  on  tQ  and  5. 

We  shall  require  several  intermediate  lemmas  before  we  are  in  a 
position  to  prove  Theorem  (8.5).  We  shall  adopt  the  following 

notation  through  the  remainder  of  this  appendix: 

Let  F and  G be  nxn  and  nxm  matrices  over  with 

sup  DF(t)n  , sup  l!G(t)il  < M 
t t 

If  for  some  t in  Z a mxn  matrix  L(t)  is  defined,  then, 

F(t)  :=  F(t)  + G(t)L(t)  . 

Also  define  for  t^  > tQ,  the  state-transition  matrix 
4.(tj,tQ)  :=  F(t^-l)  ...  F(tQ+l)F(tQ)  . 

Let  N > 0 be  a (fixed)  integer  and  let  e > 0 be  a (fixed)  real 

number. 


We  first  have  the  following  simple  result. 
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(C.2)  LEMMA.  Let  S = {xj^,X2»  X|^}  be  a given  set  of  vectors  in 

R^.  Then,  there  exists  a subset  X = {x.  , x.  . ....  x.  } of  S such 

j 1 j 2 J p ■ “ 

that  the  following  two  conditions  are  satisfied: 

(3)  The  vectors  in  X are  linearly  independent. 

(b)  Given  any  x.j  in  S - X.  we  can  write 

”1  “ *1,S  ’'J3  ’ ‘ 1 • 

PROOF.  Extend  S by  {b^.b2.  ....  b^_^}  to  form  a basis  for  R^.  Let 

B = [b]^  b2  ...  b|^_p]  and  let  S = [xj^  Xp  ...  X|^].  Define  the  full 
rank  matrix  A by 

A :=  [B  S]  . 

In  A.  select  an  nxn  submatrix  C of  largest  determinant.  Clearly.  B 
is  contained  in  C,  and  C is  invertible.  Let  C = [B  x].  and  let 
X = {x  j.x  j,  ...  X j}.  The  subset  X of  S satisfies  condition  (a) 
above.  Given  any  x^-  in  S - X,  we  can  write  (uniquely)  using  Cramer's 
rule. 


"i 


r 

E 

k=l 


det  C, 
k 

det  C 
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where  C|^  is  the  matrix  C with  k-th  column  replaced  by  x^- . Since  C|^ 
is  also  a submatrix  of  A and  since  det  C is  maximum  (among  the 
submatrices  of  A),  it  follows  that 


completing  the  proof. 


□ 


Let  X = {x2,X2>  x^}  be  an  ordered  set  of  independent 

A ^ A A 

vectors  in  R . Let  X = {x^,  X2»  x^}  be  an  (as  yet  unspecified) 

rearrangements  of  X.  Using  the  Gram-Schmidt  orthogonal i zati on 
procedure,  we  can  (uniquely)  write 


x^  = b^  + b^ 

A 

^3  “ “3,1  * ’^3,2^2  * 


V ■ “r,l»l  * “r,2'’2  * 


+ a ,b  1 + b 
r,r-l  r-1  r 


where  {b^,b2,  b^}  is  an  orthogonal  set  of  vectors.  The  vectors 

b,j  and  the  coefficients  “-j^j  very  much  a function  of  the 

A 

particular  rearrangement  X chosen.  We  now  select  a special 

A 

rearrangement  X of  X as  follows. 

A 

Chose  x^  in  X such  that  itb^B  (=  DXj^ii)  is  a maximum.  Then, 

^ A 

select  X2  in  X - {x^}  such  that  nb2H  is  a maximum.  At  the  i-th 
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step,  we  choose  x^.  in  X - {x^,X2>  such  that  itb.  n is 

a maximum. 

We  shall  now  prove  that  for  this  special  rearrangement  X, 


(C.4)  |a  I < 1 , j = 1,2,  i-i  . 

' 9J 

For  a fixed  i. 


'‘i  = " 


* “l.i-lN-l  * ^ 


Since  was  chosen  such  that  nb^n  was  a maximum,  it  follows  that 


nb. 


2 "2 
I > nx^.  n 


= z a.  iib  .n^  > a?  nb.  ii^ 

l*J  J t»i  i 

2 ^ 
Therefore,  “ij  < 1-  In  general,  for  j = 1,2,  ...,  i-1,  x.  was 

^ J 

chosen  such  that  nb.n  was  a maximum.  Hence, 

J 


(C.5) 


iib.n^  > 
J 


j-1 

E 

s=l 


2 


i 

= z 
s=j 


? ? 
nb  n > a. . . 
s i,j 


nb.n^ 

J 


2 

Therefore,  a^.  ^ ^ < 1,  proving  (C.4).  Notice  also  that  with  j = i -1 
in  (C.5) 
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2 7 7 

nb.  T n > E or.  Jib  ii^ 

’-1  s=i-l  = 

2 2? 

> a.  nb,  = nb. 

• 5 I J.  1 

Consequently,  it  follows  that 
(C.6)  ab, a > ab„a  > > ab  a 

i i r 

Let  X = {xj,X2»  X|^}  be  an  ordered  set  of  vectors  in  Define 

r(X)  = the  number  of  linearly  independent  vectors  in  X.  The  set  X 
will  be  called  wel 1 -ordered  if  and  only  if  the  following  conditions 
are  satisfied: 


a)  {X]^,X2*  Xp}  are  linearly  independent, 

b)  for  i = r+1,  ...,  k,  we  can  write 


X 


i 


r 

2 a.  . X . , wi th  la.  . I 
j=l  i.J  J ' 


< 1 


c)  on  applying  the  Gram-Schmidt  orthogonal ization  procedure  to 
•••»  Xp}  as 


Xi  = bj 


”2  ■ “2,l*’l  * *’2 
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V = “r.l  *>1  * »r,2  '’2  * •••  “r.n-l  ^-l  * ^ 


with  {bpb^,  ...,  b^}  orthogonal,  we  have 


<1 


nb,n  > iib„i  > > b n 

12  r 


We  can  summarize  the  above  discussion,  along  with  Lemma  (C.2)  by 


the 


(C.7)  PROPOSITION.  Every  ordered  set  X = {xi,x^,  ...,  x^,}  of 

vectors  in  has  a well -order  rearrangement  X = (x^.x^.  ...,  - 1 I 


For  a well-ordered  set  X,  define 


x(X)  :=  ; such  that  nb.  D -4  > n b.  , n 

1 i^N  1+1 


Clearly,  X < r < k.  We  are  now  in  a position  to  prove  the  key, 


(C.8)  PROPOSITION.  Let  t be  a (fixed)  integer  0 < t < N.  Let 

^t  ~ be  a set  of  vectors  in  R^.  Suppose 

that  there  exist  u.j(t),  i=l,  2, 


...,  k such  that 
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nu.(t)ii  < M 


(C.9) 


n<|>(N,t+l)x^  (t+1)  n < £ 

Xi(t+1)  :=  F(t)x.j(t)  + G(t)u.j(t).  Then  there  exists  an  mxn 

matrix  L(t)  with 

(C.IO)  llL(t)ii  < IL_|!L 

and  such  that 


PROOF.  Without  loss  of  generality,  assume  that  is  wel 1 -ordered. 
Let  b^- , “-j^j  be  as  in  the  definition  of  well-ordered  sets  (see 
discussion  preceding  Proposition  (C.7)).  Let  X = X(X^).  Extend 
{xj^(t ) ,X2(t ) , ...,  X;^(t),  b^^j(t),  ...  b|^(t)}  orthogonal  ly  by 


tbk+i(t),  ...,  bp(t)}  to  form  a basis  for  R'^.  Define  L(t)  by 


L(t)  : x.(t)  u.(t)  i=l,2,  ...,  X 


(C.12) 


L(t)  : b. (t)  > 0 i = X + 1,  X + 2,  ..., 


We  first  prove  (C.IO).  Notice  that 
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n E 3.b(t)ii^  = E 3?  nb.{t)n^ 
i=l  ^ i=i  1 1 


(C.13) 


^ 2 
i=l  ^ 


We  now  show  inductively  that  for  i=l,  2,  ...,  X 


(C.14)  llL(t)b.(t)ll  < 2''"^  M < 2"^  M 


For  i=l, 


nL(t)bj{t)n  = iiL(t)xj(t)o  = nu^(t)ii  < M 

and  the  assertion  is  clearly  true.  Assume  that  the  assertion  holds 
for  i=l,2,  s.  We  then  have 

a(t)bs^l(t)n  < nL(t)x  (t)n  + e |a  . | •iiL(t)b  .(t)ii 

j=l  ^ J 

By  (C.4),  l«s+i,j|  < 1 and  Hu^^^(t)ii  < M by  (C.9).  The  above 
equation  then  becomes 


llL(t)b^^^(t)ll  < M(1  + z 2'^'^)  < M 2^ 

j=l 


140 


completing  the  induction. 

Combining  (C.13)  and  (C.14)  we  see  that 


IIS  e.L(t)b.(t)«‘ 

iiL(t)n^  = sup 


u 2:  e.  b.(t)ii^ 
i=i  ^ ^ 


< sup 


( E |0  12'’  M) 
i=l  ^ 

^ 2 


.2N 


< sup 

^r**^x 


n - N 2 ( ^ lej]^ 

''  MM*'  4-1  I 


i=l 


mmN  2 

< ^n2  MM  j 


proving  (C.IO). 

We  now  prove  (C.ll).  Notice  first  that  for  i=l,2, 
F(t)x.(t)  = F(t)x.(t)  + G(t)L(t)x.{t) 


X.  (t+1) 
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Therefore  from  (C.9), 


(C.15)  H')>(N,t+l)F(t)x^.  (t)  II  < e , for  i=l,2,  \ 


Also,  for  i = X+1,  r from  (C.13), 


Xi(t)  = *"  “i 


where  Hc(t)n  < — Using  (C.3)  we  can  rewrite  the  above  equation  as 
M 


Xi(t)  = Y^^ix^(t)  + + ...  + Y.^^x^(t)  + c(t)  , 

where  |Y-jj|  < n.  Consequently,  for  i = X+1,  ...,  r, 

(C.16)  "<l>(N,t+l)F(t)x.  (t)  H < Xne  + e = (Xn+l)e 

For  i = r+1,  r+2,  ...,  k,  it  follows  from  the  definition  of  well 
ordered  sets  that  we  can  write 


’‘1“’  ' “1.J  ‘ 


Therefore,  for  i = r+1,  ...,  k. 


(C.17) 


D<J.(N,t+l)F(t)x.(t)n  < Xe  + (r-X)(Xn+l)e  < n^e 
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Summarizing  (C.15)  - (C.17)  we  see  that  for  i = 1,2,  ...,  k. 


ll(|)(N,t+l)F(t)n.  (t)n  < 


i f X=k 
otherwi se 


□ 


(C.18)  PROPOSITION.  Suppose  that  there  exist  u.j  (t),  i = 1,2,  ..., 
t = 0,1,  ...,  N-1  with 


Bu^.  (t)U  < M 
and  such  that 

Bx^ (N) B < e 

where  x.j(N)  is  defined  recursively  by  x^(0)  = e.;  (the  i-th  unit 
_vector  in  R'^,  x^(t+l)  = F(t)x.j(t)  + G(t)u.j(t).  Then,  there  exist  mxn 
matrices  L(t),  t = 0,1,  ...,  N-1  and  a real  number  A with 

BL(t)B  < A 
and  such  that 


BF(N-l)F(N-2)  ...  F(1)F(0)B  < . 
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PROOF.  We  first  set  up  some  notation.  For  a set  of  vectors 
X = {xj^,X2,  ...,  X|^}  in  R^,  let  X = {x^,X2»  ....  X|^}  be  any  well- 
ordered  rearrangement  of  X.  Recall  that 

A A 

X = x(X)  :=  ; such  that  nb.  n > -^  > iib.  ^ « , 

1 [v|IN  1+1 

and  that  k = k(X)  = the  number  of  vectors  in  X.  Define  a set  of 
i ntegers 

A A A 

1^  •“  {j  • ^Xj^,X2»  •••»  * 

For  t = 0,1,2,  ...,  N-1  define  X^  and  recursively  by 
Xo  = n 

Xq  - {Xj ( 0) ,X2 ( 0) , ...,  x^(0)}  , 

\+l  ""  ^X-j(t+l)  : i£  } 

From  the  above  definitions  it  is  clear  that 


(C.19)  n = Xq  = k 


1 


> X^=  k2  > X2 


\ = 


't+1 


> X 


t+1 


> X 


N-1 


> 0 


• • • 
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Define  f^,  t = 0,  1,  N-1  by 


f 


t 


1 


4 


n 


i f X, 


otherwi se 


It  is  clear  from  (C.19)  that  (and  this  is  the  key  step) 

(C.20)  ^N-1^N-2  •“  '*'0  ^ • 

We  now  apply  Proposition  (C.8)  to  to  conclude  that  there  exists 

n ^ 

an  mxn  matrix  L(N-l)  with  HL(N-l)n  < CM^  where  C = and  such 

e 

that 


»F(N-l)x.  (N-l)ll  < efj^_^  , X.  (N-1)  e . 

Notice  that 

IIF(N-1)D  < HF(N-1)D  + nG(N-l)L(N-l)D 
< M + M(Cm'^)  < 

We  can  again  apply  Proposition  (C.8)  to  X|^_2  (with  M being  replaced 
by  CM  ^ and  e being  replaced  by  £f|\j_2)  to  conclude  that  there  exists 
an  mxn  matrix  L(N-2)  with 
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IIL(N-2)H  < (CM)^ 
and  such  that 

llF(N-l)F(N-2)  x.(N-2)»  < . x.(N-2)sX^_2  . 

Repeating  this  argument  N times,  we  conclude  that  there  exist 
matrices  L(0),  L(l),  L(N-l)  with 

nL(t)n  < (CM)'^  =:  A 

and  such  that  for  i = 1,2,  ...,  n, 

IIF(N-l)F(N-2)  ...  F(0)e.ll  < • . . f „ . 

/ 

The  above  equation,  along  with  (C.20)  gives  us 
nF(N-l)F(N-2)  ...  F(0)n  < n®e  , 
completing  the  proof. 


We  are  finally  in  a position  to  prove  Theorem  (8.5)  which 


mxn 


□ 


is 


restated  below  for  convenience. 
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(8.5)  THEOREM.  Let  Z - (F,G,H)  be  a linear  time-varying  system  over 
^”(Z).  Then  I is  asycontrol Table  if  and  only  if  Z is  stabilizable. 

PROOF.  We  have  already  (see  discussion  preceding  Theorem  (8.5)) 
shown  that  stabilizability  implies  asycontrol 1 abi 1 ity . We  now  prove 
the  converse: 

Suppose  E is  asycontrol lable.  Choose  in  the  Definition  (C.l)  of 

asycontrollability,  e = — Then,  by  Proposition  (C.18),  there 

2n'^ 

exist  matrices  L(t),  t = 0,  1,  ...,  N-1  such 
nL(t)ii  < 2An^ 


HF(N-l)F(N-2)  ...  F(l)F(0)n  < ^ , 

where  F(t)  = F(t)  + G(t)L(t).  Recall  that  here  N is  the  number  of 
steps  required  to  drive  all  states  to  energy  < e as  in  Definition 
(C.l)  of  asycontrollability.  Since  all  our  arguments  are  independent 
of  initial  time,  we  can  (again  by  Proposition  (C.18))  find  matrices 
L(t)  for  all  t in  Z such  that  HL(t)n  < 2An®  and 

llF(kN-l)F(kN-2)  ...  F(kN-N+l)F(kN-N)  n < -^ 


for  all  k in  Z.  Thus  the  matrix  L(t)  view  as  being  over  A is  bounded 


(i.e.  in  a“(Z)),  and  (zI-F-GL)"^ 


Thi  s 


□ 


completes  the  proof. 


is  a stable  power  series. 
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